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SYNOPSIS 


Thesis entitled ^ELECTRONIC STRUCTURE OF DISORDERED 
ALLOYS' submitted by VIJAI KUMAR SRIVASTAVA to the 
Department of Physics/ Indian Institute of Technology 
Kanpur, in partial fulfilment of the requirements of the 
Ph.b. degree. 


Coherent-potential approximation (CPA) , during the past 

two decades, has been one of the main theoretical developments 

in the area of study of electronic st3ructure of disordered 

alloys. The simplest case for which most of the CPA calcula- 

tions have been done, is that of a random binairy alloy 

A B, where c is the concentration of A in B. Here the 
c 1-c 

site energies are taken to be random (diagonal disorder) and 
the overlap- integrals are assumed to be non— random, i.e* I 

there is no off-diagonal disorder. Also, it is further ' 

assumed that correlated scattering from two or more sites is 
ignored. This is known as single site coherent-potential 
approximation (ICPA) . There are several equivalent approaches 
to ICPA (reviewed in the chapter II of the thesis) all of 

which lead to the same results. | 

[ 

The need to go beyond ICPA has been felt for a long time. 
Several suggestions have been put forward about the ways of 
generalising the ICPA and ways of incorporating off-diagonal 
disorder and short-range order (SRO), Calculations based on 
these suggestions suffer from two main drawbacks. Firstly, 
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for any exactly solvable realistic model, which involves 
more than one orbital per site, the corresponding algebraic 
eguations become tediously unwieldy, particularly on a 
three-dimensional lattice. Secondly, any tractable or further 
simplified approximation always leads to the unphysical 
averaged Green's function in the strong scattering regime. 

The essential Herglotz property of the Green's function is 
violated. 

The correct Herglotz property of Green's function is 
preserved in the augmented- space formalism introduced by 
Mookerjee (1973). This formalism is discussed in detail 
in chapter III. The present work which is based on augmented- 
space formalism offers an unambiguous, tractable, self- 
consistent cluster generalisation (CCPA) of ICPA, which 
retains the Herglotz property at all energies for all ranges 
of disorder. Off-diagonal disorder and short-range order 
have also been incorporated. The general fojrmalism of the 
CCPA is also discussed in chapter III. A graphical method 
coupled with the recursion technique of Haydock (1972) enables 
us to work on a realistic three-dimensional lattice. 

The most important contribution of this work is the 
generation of a self-consistent self-energy. This reproduces 
the density of states near the band edges more accurately 
than that produced by non, self-consistent cluster embedding 
approaches tried earlier (Choudhary, V. 1981) • The embedding 
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in a non self-consistent medium produces spurious structure 
near the band edges, which is completely absent' here. 

The results of a self-consistent cluster calculation 
of a single orbital model on the diamond and b.c.c. lattices 
have been reported. The most instructive feature of the 
results is that here the band width is larger than that in 
ICPA. Clustering features appear in the minority (impurity) 
band and effects of off-diagonal disorder and SRO are easily 
identified. 

In the chapter IV, it is shown that how our cluster-CPA 
formalism can be applied for the calculation of electronic 
density of states of valence bands of ternary, random III-V 
semiconducting alloys e.g. Ga Tu, _ As where 0 < C < 1. 

C X *“C 

3 

In these alloys each atomic site provides four sp hybridised 
orbitals for bonding. These materials form tetrahedrally 
cordinated, covalently bonded solids in a zinc-blende 
structiire. In these alloys one f.c.c. sublattice of the 
common ion is non-random while that of the other ion is 
randcmj^ly occupied by the constituent atoms. 

Chen and Sher (1978) have applied single- site, 
self-consistent CPA to the III-V semiconducting alloys. 

They used a bond orbital model (BOM) and parameterized 
Hamiltonian fitted to experimental data. This approach 
is extremely instructive but completely ignores the effect 



of clusters and cannot adequately deal with off-diagonal 
disorder. They have treated the diagonal-disorder in- 1 CPA 
but off-diagonal disorder in virtual-crystal approximation 
(VGA) . This is not a consistent procedure. Here an 
attempt has been made to improve upon this drawback. The 
coherent-potential medium is made self-consistent and off- 
diagonal disorder is taken into account. A calculation of 
this type has been done for Ga„ ^In„ „ As. The interaction 
between adjacent bonds with an anion in common has been 
taken into account exactly. The result is that there is an 
improvement in the band-edges and structure near them, 
particularly in the p- like part of the density of states. 

In the concluding chapter V, the possible extensions 
and future usefulness of the cluster-CPA have been 


discussed 



CHAPTER I 


INTRODUCTION 

1.1 Disordered Solids: Classifica-bion and Problems 

Ever since the enunciation of the Bloch^ s theorem in 
1928 and the emergence of the band picture, considerable 
successful effort has gone into the understanding of 
electronic properties of crystalline, periodic systems. 
Concerted and extensive effort began into the study of 
disordered systems only in the late sixties. However the 
following decade has seen extensive development in what was 
then termed. a ^frontier area* of research. Experimental 
as well as theoretical understanding has undergone such a 
seachange, that disordered systems are no longer considered 
' exotic * . 

Disordered solids encompass a rather extensive class. 
Broadly, they can be grouped into following classes: 

(1) Compositionally disordered -solids as in the case of , 
random alloys both substitutional (CuZn, Ga^ In^_^As etc. 
and interstitial (PdH) . 

(2) Structurally disordered solids as in the case of 
amorphous materials like a-Si, a-Se, a~GaAs and 
glasses etc. 

(3) Magnetic alloys with compositional disorder (spin- 
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glasses/ AuFe, CuMri/ Fe^!sli^_^/ Fe-Ni-Mn) and structural 
disorder (e.g. amorphous magnets like metal glasses ). 

(4) Amorphous alloys e.g. a-Ga^In^_^As/ a-Si(H) . 

(5) Polymers. 

The kind of problems one encounters when one turns to 
random systems are also multifarious . For example/ 

(1) The understanding of various disorder linked transitions/ 
magnetic transitions/ order-disorder transitions. 

(2) The theoretical understanding of various electronic 
properties like density of states, optical properties, 
photo- emission, conductivity, dielectric response and 
electron related magnetic properties. 

( 3) The understanding and description of the nature of 
disorder- driven localization-transition of the 
excitations in random systems. 

(4) The understanding of other elementary excitations 

in random systems such as phonons, magnons, polarons etc. 

In the follov/ing work we shall concentrate on compo- 
sitionally disordered crystalline alloys. Further we shall 
address ourselves exclusively to the londerstanding of density 
of electronic states. Although the various techniques we 
shall introduce, may with a little modification also be 
applied to other classes of systems and certainly to other 
electronic, phononic, magnonic etc. properties. We shall 
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not touch upon phase-transitions of any kind nor investigate 
the ^Localization^ problem. 


1 .2 Theoretical Approach 

The structure and properties of various crystalline 
materials can be very well understood by the application of 
Bloch's theorem. The idea of bands and band gaps came out 
of this- But/ the long-range order that Bloch's theorem 
assumes/ is not essential for the existence of bands and 
gaps between them. The well defined bonding and antibonding 
energy levels in isolated atoms spread into bands as the 
atom themselves bond to form the solid state. This has 
been the view point of chemists. In fact all the band 
structure calculations of density of states take account 
only of interactions between near neighbours. Long-range 
order is not necessary. The local bonding behaviour is 
sufficient to determine the energy bands and properties 
arising from these. This fact suggests that a theory based 
on local- chemistry of bonding may be attempted. 


The usual starting point is Schrodinger equation: 


H If) = E ,1, 
— -n n^n 


( 1 . 1 ) 


However the eigenfunctions are inappropriate 

for the discussion of large systems with consequent high 

degeneracy because any physically measurable quantity 

23 

involves summation over N 'v 10 electrons with a loss of 
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physical insight and anything related to the eigenfunction 
Also the eigenstates are unstable under small pertur- 
bations which are insignificant beyond some distance from 
perturbation. The k-space, which relates to the individual 
eigenstates (which are inappropriate for our purpose) / 

needs to be ''thrown out'' (Heine, Solid state physics, 
vol. 35, p« 5). This can be achieved by taking recourse to 
what is called 'local density' of states defined as: 

n(E,r) = n(E) lif^gCr)!^ (1»2) 

where n(E) is the total density of states of the system, 

2 

that is modulated by | i for a typical state of energy E« 
Integrated upto Fermi-energy Ep, this gives the total 
electron density. In an alloy or compound it tells which 
bands are more or less weighted on which atoms. 

Alternatively, for each spin direction one can write: 

n(E) = I 6(e-E^) (1.3) 

n 

n(E,r) = I it')^(r)|^ 6(E-E^) (1.4) 

n 

Here the summation is over all the eigenstates of H. 

Equation (1.4) defines the local density of states 
interms of extended eigenstate 4'j^ of the whole system and 
its determination via the determination of and E^ 

(i.e. band structure) helps us in no better way. But 
fortunately n(E,r) is related to the imaginary part of 



of the Green's function by the relation (see appendix B), 

n(E,f^ = - 1. im G(r, r/E; (1.5 ) 

Jl 

and the Green' s function can be obtained directly without the 
knov/ledge of ’4' , The Green's function fomaalism provides 

"" , cJI 

one of the most powerful methods, since ef the electronic 

properties and their response to perturbations caused by 
externally ajpplied forces are related to the single or many- 
particle Green' s functions. 

We shall adopt a tight-binding approach as most suitable 
for a local environment view-point. We may express eigen- . 
states as a linear-combination of atomic orbitals (LCx\0)s 




= I 




(1.6) 


Here 2. is the site index for some atomic site and O- indicates 
the type of orbital. The orbitals 6 (r) are localised 

along the atomic site i (but this has no bearing on the fact 
that are extended). In view of equation (1,6) H can 

be written as a matrix with elements <a 2, i H i a' £' > ( in Dirac 
notation) representing the hop^ping integrals of the electron 
from orbital al to a' £' for a'£' ^ a£ and the energy of 

the orbital a £ for a' £' = a £ , 


The principal difficulty with the old LCAO methods vjas 
the non-orthogonality of basis i.e. 


<a £ia' £'> 


,aa' 

££'' 


7 ^ 5 


££' 
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This leads to computational difficulties. There have been 

various attempts at methods which '’ignore' the non- orthogonality 

occx ^ 

and assume S = g § and take into account the 

aa' 

effect of non-orthogonality in different, not very satisfactory, 
ways. Other efforts change the basis from |a£> tc U|a£> 
such that U’^’SU = This new basis is orthogonal and may be 

used directly. (Anderson, .1969) . 

An old approach resuscitated in recent years has provided 
a real break-through in making the LCAO method feasible and 
eventually as good as any k- space methods-this is the 
''chem.ical pseudo-}potential' ' method. (Bullett, Heine, 

Solid State Physics, Vol. 35). The basic idea of the method 
is the use of a new pseudo-basis 1§ and a new effective/ 

v-L 36 / 

~ —1 

Hamiltonian H = S H, both of which may be obtained from the 
atomic wavefunctions and potentials (Hermann and Skillmann, 

1963) . Here H is non- symmetric but that in no way creates 
difficulties in the procedure v/e shall subsequently follow. 

In some of the older books LCAO is described as ultimately 
Vinfeasible because to use it effectively the computational 
effort becomes intractable. With the introduction of 

chemical pseudo-potential method, various recursion techniques 
(Haydock, Heine, Kelly, Solid State Physics Vol. 35) and fast 
computers, this drawback no longer exists. LCAO or modifi- 
cation thereof, has had a renaissance (Bullett, Solid State 
Physics Vol. 35. p. 129). 
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The resurgence of the tight-binding/ local environment 
methods is a great boon to us — as in disordered systems the 
lack of periodicity invalidates all k-space approaches. 

Assuming now that basis-orbitals are orthogonal, it can 
be easily shown that the local density of states is given by 

= " I <a£iG(E)|a£> (1.7) 

where G(E) = lim ^ (see appendix B) 

Z-*E 

Here the Green's function G(E) is a matrix which can be trans- 
formed according to different basis functions. The important 
point to note is that for determination of n^j(E) we need 
only one matrix element of G(E) • The same is tme for 
equation (1.5); we need to know only r = r' part of GiE) 
which simplifies the problem enormously. We are able to 
express measurable quantities in terms of some appropriate 
small parts of G(E) that can be solved for and computed 
seperately from all the unwanted remainder of G(E) . 

Local orbitals viz. 'the atomic site orbitals or the 
bond-orbitals in case, of covalently bonded atoms and their 
interaction between near neighbours only, have been used as a 
basis for small as well as large systems to explain the 
observed phenomena. In alloys also, affinities are known to 
affect essentially the near neighbour configurations. This 
is a pointer to a valid local— theory and is welcome because of 
its implication that considering configurations of local 
environments i.e. small clusters, may be sufficient. 
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In the absence of disorder# once the Hamiltonian model 
is set up# the solution of cjuantum-mechanical problem is 
essentially the solution of Schrodinger equation. However, 
with the introduction of disorder we face a new problem i.e. 
the lack of sufficient inforrriation about the sample. For 
example, in a random, binaary alloy, the configuration of 
individual atoms on the lattice sites is different for 
different samples and an experimentalist at the outset has 
no knowledge of the particular configuration of his sample. 

Nor is such a knowledge of much use. Evidently, every 
different sample of the alloy will give experimental information 
having different micro- structure according to the particular 
configuration it has. The experimentalist, on the other hand, 
is primarily interested in the overall statistical trends 
of his results which he interprets as the physical 'properties'^ 
of the random alloy. A theorist too can generate a number 
of results depending upon the particular configuration he 
chooses for his random parameters. This again is not very 
useful. One could try to theoretically derive properties 
of the most probable configuration and propose that these 
properties are highly likely. But again the number of 
different configurations as probable as the most likely and 
those of comparable likelihoods of occurrence, are intractably 
numerous. Therefore, v/e have to evolve a statistical 
discription of the system. 
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The most basic of these descriptions is the averaging 
over all possible configurations. It was pointed out by 
Anderson in 1972 that the type of statistical description 
essentially depends on what v/e are trying to describe. He 
warned that only experimentally observable physical quantities 
related to the response of a physical system to a probe must 
be averaged. This vjraming was given to explain source of 
controversies that arose over non-physical results obtained 
because non-physical quantities like self-energy, the wave- 
functions etc. were averaged. The Greenes function corres- 
ponding to any spectrum like that of the electrons, phonons 
etc., can be configurationally averaged since its imaginary 
part is directly proportional to the density of states, a 
measurable quantity. Other physical quantities may be 
derived from it and compared x-/ith experiments. 'This is an 
additional pointer towards a local theory. 

In fact, the study of configurationally averaged single- 
particle Green's function forms the substance of our thesis. 

1.4 The Model 

The full many-body electronic problem in the disordered 
system is beyond the scope of present work. We shall consider 
an equivalent one-electron problem. The one-electron 

potentials will be chosen in such a way as to incpfrporate 

/ 

the maximum possible many-body effects. 
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Firstly, the Bom-Oppenheimer approximation helps us 
to separate the nuclear and electronic parts of the vjave— 
function describing the many-particle system consisting of 
nuclei and electrons in the alloy. The Schrodinger equation 
is still a m.any electron v^ave equation describing the motion 
of electrons in a field of nuclei. The Hartree-Fock scheme 
then yields a one- electron equation for stationary states 
in the presence of nuclei and an average potential due to 
all other electrons. The effects of exchange and correlation 
can also be taken into account by Slater^ s method. 

In case of a random, binary alloy of the type , 

the Hamiltonian in a Hilbert space can be written ass 


H 


.1 


e . 
xn 


P. 


—in 


+ . 1 . 


I 


y V. . T. 

^ i5i,jm -in,jm 


( 1 . 8 ) 


Here P and T are projection and transfer operators, 
i and j are nearest neighbour site indices and n and m 
are band or orbital indices. Here both 

£V. . } can be random variables. The disorder, due to 

in,jm 

randomness of e. , the atomic potential at the site i, 

in ^ 

is called diagonal disorder whereas the disorder due to 

randomness of V . . , the overlap- term, is called off- 

i.n,jm 

diagonal disorder. Here e^^ can take two possible values 

A 3 

or depending on whether the i-th site is occupied 

by A-t^rpe or B— type constituent. The dependence of e^^, 
on the configuration around the i-th site is to be ignored. 
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The occupancies of each site are assumed to be statistically 

independent. Short-range order due to chemical affinities 

is not taken into account. The probabilities of finding an 

A-type or B-type constituent at any site are then directly 

proportional to their concentrations. Here V. . can 

in^jm 

take three possible values , V^~ or 

^ nm nm nm 

The main idea is to get a configurationally averaged 
Green ''s function i.e. <G(E)>, incorporating the effect of 
clusters. In viex^ of this, we have first reviewed in 
Chaxoter II, the older methods (i.e. single-site coherent- 
potential approximation methods) of calculating <G(E)>. We 
have pointed out the drawbacks and discussed possible lines 
of generalisations. In Chapter III, we have developed a 
self-consistent cluster-coherant-potential approximiation 
(CCPA) based on augmented space formalism of Mookerjae (1973). 
This formalism (CCPA) has been applied to three-dimensional 
diamond and b.c.c. lattices (Chapter III) choosing arbitrary 
values of atomdc potentials and overlap-integrals. The 
results are quite encouraging. The CCPA formalism has been 

applied to III-V semiconductors and their alloys using bond- 
orbital model (Chapter IV) . These semiconductors and their 
alloys are covalently bonded with a four- fold co-ordination 
and the alloys fall in the category of substitutionally 
disordered solids. 

In the concluding Chapter V, we have discussed the main 
achievements of the CCPA fomalism and its various possible 
applications. 



CHAPTER II 


OLDER FORMALISMS AND THEIR DRAWBACKS 

2.1 Introduction : 

In this chapter we shall review the existing methods 
for calculation of the electronic density of states in disordered 
solids. The aim will be to examine these methods in detail with 
a view to point out the approximations involved in them and 
particularly the shortcomings that these approximations entail 
in the results. Such an examination will indicate a way of 
generalising the methods so as to overcome the shortcomings and 
formulate a method that is feasible and yields physically valid 
results in a realistic situation. 

There are several approaches available to us of increasing 
degree of sophistication such as Rigid-Band approximation/ Virtual- 
Crystal approximation/ Minimum-Polarity approximation. Virtual 
Bound State approximation and the Coherent— Potential approximation. 

Except for the single site coherent-potential 
approximation (ICPA) and virtual-crystal approximation (VGA) /all 
other approaches mentioned above have very limited success and 
are suited for some specific system and fail in others* We shall 
review these two in turn* 



13 


2*2 Virtual-Crystal Approximation : 

In this raodel the actual random alloy potential is 
replaced by an average, periodic potential 

V = cV^ + (1-c) Vg 

x^hers and are the potentials of the constituents A and B 
present in the atomic concentrations c and (1-c) respectively. 

Tl'ie VCA does not reveal any features characteristic of disorder 
other than the crude averages. In case the potentials of 
constituents are very similar, the VCA is not a bad approximation. 
Because the effective potential is periodic, the standard methods 
available to carystalline systems are valid here. The VCA often 
serves as the starting point of more sophisticated iterative 
self-consistent approximations. 

2.3 Coherent Potential Approximation (CPA) : 

The single site CPA had remained till late, the main 
theoretical development in the study of electronic stxncture of 
random systems. It maintains the analytic features of the exact 
Greenes fiinction and interpolates correctly between several 
individual limiting cases e.g. virtual-crystal, atomic and dilute 
limits , in the case of random, binary alloys. 

The literature on ICPA can be classified into three 
categories : (1) Diagrammatic approaches using propagator or 
locator formalisms (2) The ef fective-mediim approach (3) Multiple- 
scattering approaches. All different approaches yield the sanie 



final result x^ithin this approximation. A detailed account of 
all the three approaches is given here. 

2.3.1 The Diaqramrrtatic Approach ; 

This approach was introduced by Yoneza\'i?a and 
Matsubara (1966)/ Yonezav/a (1968) and Leath (1968,70) using 
propagator formalism. VJe shall illustrate the approach through 
tight-binding h.nderson model (1958) including diagonal-disorder 
only. The Hamiltonian can be written as s 


H = Hq + V 


( 2 . 1 ) 


Here is the non-random translational ly invariant part ^ jj ' 


while V = Iv. = I e.P. is the random part of H. ' e. is 

- . -1 . 1—1 - 1 




random variable X'^hich can take values e^^ or eg where e^ and eg 


are the atomic potentials of the constituents A and B respectively 


and T^ . are projection and transfer operators. The unperturbed 

—I 

propagator is G = (ZI - H^) . The Dyson equation relating 


■ o ■ - -o 

the perturbed propagator G to G is 


G = Gq + Go V G 


= So + So Y So + Y Y So + •• 


( 2 . 2 ) 


On configuration averaging either side \-j±th respect to variables 


{ e^ } we have 


<G^ = G_ + Z G ^ G_ + Z G *7 "'^4^ G +. . . 

— ■ — O O ““1 —O . 'O — i—o — J -"O' 

1 1/3 


(2.3) 


Here <G> represents configuration average of G. Equation (2.3) 
forms the basis of diagrammatic approach. For the generation 
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of diagrams# for every propagator a horizontal line is 
assigned. For every <G> we assign a double horizontal line 
and for each within the averaging# we assign a vertical 
dashed' scattering-line' # together with an 'interaction vertex' 
labelled i (see fig. 2.1). Terms like G_ <v. G_ v. > G_ are 
represented by 

x^ ■ ' , - 

/ - X 


The diagrammatic representation of equation (2*3) is shown in 
fig. 2.2 - i;. 

We call a diagram 'irreducible' if it cannot be broken 

~ , >c 

into two by mpturing along a propagator. An irreducible 
diagram is called an 'skeleton' diagram if it does not consist 
of units enclosed within one another# so that the enclosed 
units could be considered as insertions on the internal 
propagators* In fig. 2.2 diagrams (a)# (b)# '(d)# (h) and (i) 
are irreducible of which (h) is not a skeleton diagram 


In order to sum the contributions of various diagrams 
a set of rules# similar to Feynman rules# is evolved’. We shall - 

‘"v ' ' ' ! 

take the specific case of a -random, binary alloy for illustration.! 
The ecruation (2.3) is resnritten as follows : 


<S> = J So -'YP So + So ^Yi So YP So 

+ y G ^ V ■ G^ V ■ ^ G_ +■ ... 

u —Q — 3_ ■•.'■O '-*1 . —O'. 


(2.4) 
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Now we introduce the following contribution irules ; 


(i) 


With each horizontal oropagator we associate a factor 
= (ZI - where 


H = I e„ P. + I V. . T. . 
-o ^ B -1 -ij-ij 


With each scattering line we associate a factor W and with 
each interaction vertex a factor c, the concentration of A 
cornoonent. These factors arise from 


<v.> = <(e--e„) P.> = c W P., here W = (e, - e„) 

— X liD— 1 An 

is the difference in atomic potentials of the two 
constituents. 

(ii) All internal vertices are summed over without restrictions. 
The latter rule necessitates multiple occupancy corrections. 
To understand this/ let us use the two rules on the 
diagrammatic expansion of fig. 2. 2. We get 

<G> = G^ + CW y G^ P. G„ + c^ T G^ P. G^ P , 

"“O " ->»0 ***X *0 ^ "“O ""X ""C •• 1 *—0 

i i j 

Gq Pj. Gq Pj_ Gq + (2.5) 

i 

Now using the independence of various e^ we have 

G^ Vj> = Gq <yj> if i j, and the 

above equation may be v?ritten as 

<G> = + I <Y.> + I G„ <v. G„ V > Q 

1 IJ 

, i' '' 


( 2 . 6 ) 
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Comparison of equations (2.4) and (2.6) reveals the 
fact that the multiple -scattering occupancy corrections arise 
because of the sumrriation restrictions in equation (2.4). We 
may rewrite equation (2.4) removing these restrictions# to get 


< 


= 9o + I So So .1 So So 2c 

1 ij 


+ I tGo 2o =o - 9o 2oi+--- <2.7) 


The subtracted term corresponds to the contribution of 
a diagram as shown in fig. 2.3a. This diagram is obtained from 
a diagram shown in fig. 2.3b by breaking away the scattering 
lines from the interaction vertex. If we look at a higher 
order term 


I ^So ^^i So S± So 2i> 2o 


32o <Yi> 


2o <Yi 


V. > 

-o -1 


G 

-o 


- So <2i> So 


<v.> 


So So> 


the subtracted terms again correspond to contributions from 
diagrams formed by breaking away the scattering lines from the 
interaction vertex in every possible way. 

Now the diagrams may be summed according to the rules : 

(i) All skeleton diagrams are drawn and from these skeleton 
diagrams derived diagrams are generated by breaking away 
the scattering lines from the interaction vertices in 
every possible way and giving diagrams already explicitly 


included 



(ii) With every rightmost and internal propagator a contributio 
G is associated and with leftmost propagator a contri- 
bution G^. For a binary alloy, W is associated with 
eve 2 ry scattering line and a factor c v/ith every intera- 
ction vertex. An overall negative sign is associated to 
every derived diagram. Pig, 2.4 illustrates the diagrams 
generated according to new rules. 

We have now obtained a Dyson equation connecting G to 
G given by 

2 = So ^o 2- 2 (2.8) 

This is not a trivial result,- Although the Dyson 
equation exists for the unaveraged G, a priori it cannot be 
said to have existed after averaging^ This is because 
different 'self-energy' parts of a diagram for a perturbation 
expansion for ^ connected only by the propagator cannot 
contribute independent matrix factors on averaging, ' The Dyson 
equation has been derived by laborious resummation and 
regrouping of terms. 

We are now in a position to make approximations. The 
single-site GPA corresponds to choosing only those skeleton 
diagrams which involve just one site. Diagrams such as shown 
in fig, 2^,5' involving pairs and triplets etc, are rejected. 
However for self-consistency diagrams such as shown in fig, 2 . 5 a 
are included. 

The diagrams are now tabxilated as shown in fig. 2,6, 

The first column includes all skeleton diagrams, each row having 
a diagram of one order, - The derived diagrams from it occupy the 



(b) 

Fig. 2-4 

Construction of derived diagrams 



(a) (b) 


Fig. 2-5 


Ill- 



Fig. 2- 6 Diagrammatic calculation of self-energy by Leath's method 
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same row. The second colmnn contains all irreducible derived 
diagrams and subsequent columns contain derived reducible 
diagrams that break up into two, three etc. irreducible parts. 

We may nov/ add up column v^ise.' 

In ICPA the self-energy S is site diagonal i.e. S = S x. 
The sum of the first colimn is 


(bare) = cW + c W G W + c W(G V7) + 


. s = <=>•■ 


11 


1-GW 

Now S (g) = S (bare) - {S (i) - S (g)} - S (i) g 2 (t) - 


where t = G/(l- 3(i) G) takes into account internal insertions 
in the derived diagrams. The self energy is a function of G. 

A little algebra immediately gives 

E(g) = cW/{l-(V7-2) G} (2.9) 


Equations (2.8) and (2.9) are sufficient for the determination 
of S and G. S being site diagonal we have from equation (2.8) 

G = Equation (2.9) is solved by iteration techniques, 

starting from, say, the VGA (Brouers et al 1973, Ducastelle 1972). 

This process for diagram summation is quite general. 

Tor glassy disorders we have different contributions for the inte- 
raction vertices and scattering lines. 

In equation (2.1) the alternative choice of V as the 
unperturbed Hamiltonian leads to the locator formalism, where the 
perturbation expansion is in terms of localised states 
(Matsubara and Toyozawa 1962, Matsubara and Kaneyoshi 1966, 


man 1 969 ) . 
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In the propagator formalism v^e start off with 3loch-like 
extended electronic states in a periodic system. We then intrOvduce 
the diagonal disorder as a perturbation and study the effect on 
electronic states. In locator formalism we begin with 
atbmistically localised electronic states/ introduce the overlap 
term and study the resulting change. 

For the single— site CPA both the propagator and locator 
formalisms give identical results. 

2.3.2 Mean-Field Approach ; 

The idea of mean-field or effective meditam (Soven 1967) 

gives perhaps a simpler physical insight into the CPA. Here vre 

define an effective Hamiltonian H __ such that the Greene's 

-eff 

function corresponding to it is equal to the configurationally 
averaged Greenes function of the actual Hamiltonian H under 
consideration . 

Let our actual Ham.iltcnian in the tight-binding basis be 


H = I e, + I V T 

i 1 1 


( 2 . 10 ) 


where P^ and T^^ are projection and transfer operators respectively 
Here both e^/ the site energies and the nearest neighbour 

overlaps may be random. 


Within a single— site approximation it is not very 
profitable to include effects of off-diagonal disorder. 
Conceptually random / involve two sites. Intuition expects 
the lowest order approximation to be a pair-CPA ( 2CPA) . 
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In the presence of off-diagonal disorder, the ICPA density 
of states is correct only to three raoraents rather than eight 
(as in the case of diagonal disorder only) v/hereas a consistent 
2 CPA is correct up to fifteen morients, whether we have pure 
diagonal disorder, pure off-diagonal disorder or both. 


in ICPA, V. .'s are supoosed to be non- random and we 
xj 

set V. . = V, where V is some miean potential.- 
Now we define our effective Hem.iltonian 


H 


eff 


= I S (E) P + I S . (ir.-?. |,E) T 
i - ij ^ 


ij 


( 2 , 11 ) 


Again in ICPA the self-energy is sam^e for all the sites 

i.e. E . ( E) = S (E) and S. . = V. So we write finally 


H 


eff 


I 2q(E) P^ + V I T 


IJ 


ij 


( 2 . 12 ) 


The idea is to find out those S^(E) which yield 


= <G> 


(2.13) 


<G> indicates configuration average of G 


o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

c 

o 


o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 


Fig. 2.7 
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Nov-/ in order to determine H self-consistently for 

—erf 

ICPA, we embed an exact ootential e. at the site i within 
the effective medivim in such a way that we do not produce any 
extra scattering on the average (Fig. 2.7). In mathematical 
terms we have 


<G^ > 


G 


eff 


where G*-^^ = 


and + Z Ce,-E„)P. 


H 4- h 

-eff^ - 


(2.14) 


Now we may write 


,(i) 


^ 2i~Sef 


Seff 


Let us define an operator T such that 


= Seff + Seff S Seff 


(2-15) 


then 


T = (I - hSeff^" - ^ ? 


£■ P- 
—1 


1~ £-• o_ j-j: 


•"i-efi 


= I h 


Here 


^'i = ^®i " * 


From equation (2.15) vie can write 

= £eff + Seff 5eff • 


Now in order to satisfy equation (2*14) we write 


<t.> = 0 


or 


e . 

i 


S^(E) 


> = 0 


l-[e, - E^(E)]g(Z-S ) 


(2.16) 

(2.17) 
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which is identical to equation (2.9). Sc we conclude that 
in single- site CPA mean- field or effective-medium approach 
is equivalent to the diagammatic approach. 

2.3.3 Multiple- Scat t er in O' Approach 

The calculation of <G> is carried cut easily within 
the framework of multiple-scattering theory. We express 
single-particle Hamiltonian of equation (2.1C) as 

H “ £"3^ i")! (2.13) 

where is some suitably chosen reference Hamiltonian. 

The difference W is the perturbation causing scattering. 

The basic requirement of the multiple scattering theory 
(Anderson and Mcmillan 1967) is that the perturbation term 
W can be decomposed as a simn of contributions from each 
site i.e. 

W = H — H = y w P / v/ =3 — S(-S) (2.19) 

...I. ““a. “ ii 

n 

As defined earlier 

T = W + W T (2,20) 

This relates G = (ZI - H ) ^ to G by 
— a. “"Cl 

G = G -f G T G ( z • 2 i ) 

““''a — «<— cl 

In order to express T interms of single-site operators 
we use the relation (2.19) to obtain 

T = I w^(I + T) = I t^ (say) 
n ^ n 


( 2 . 22 ) 



27 


which gives 


1 = I “nEn + ^ Sa ^ "ra £m +‘ ' • 

n n m 


(2.25) 


equations (2.22) and (2.23) yield 


t =w P (I + Gt + I G t 
n -n —n —an , —a m 

m;en 

m^n 


= T (I + G y ^ ^ 

-n -a 4/ t ) 

m/^n m 


(2 • ^ 4 ) 


—1 

where = (^“W^^Pj^G^) n atomic t-matrix 

representing complete scattering from, n-th site. 

Iteration of equation (2.24) leads to the usual 
multiple scattering expansion 


t=T+TG yi+XG IxG It+. 
-n -n -n -a -m -n -a -m -a -p 


-n — n —a 


'm/^n 


iry^n 




which reveals the fact that t^ is actually dependent on all 
sites. Prom equation (2.21) on averaging both sides, we get: 


<G> = G + G <T> G 

— • — 3 . “*“3 ““ 3 . 


(2.25) 


The single site CPA can now be introduced by averaging 


(2.24) , 


<t > = <X„(I + G^ I t )> 


— n 


-n 


•a . — m 

m/m 


= <X > {1 + I <t >) 

“■ iTV^n 

+ <(T - <T >) G I it - <t>)> 

-n nv^n 


(2.26) 
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The first term in equation (2.26) may be thought of as 
representing the scattering of an isolated site i iramersed 
in an average medium, ^^rhile the second term describes 
fluctuations around this average behaviour. In single~site 
approximation we neglect the second term. The electronic 
properties of the alloy are then determined from the closed 
set of equations? 

<t > = <1 > (I + G I <t >) (2.27) 

-“11 — — a y — m 

m/n 

The self-consistent solution of follows from <T> = 0 

“**a 

which imolies <G> = G . The ICPA eauation for determining 

— —a 

self- energy nov; becomes 

<T > = 0 

— n 

(e - S (E)) 

or < — -bi > = 0 (2.28) 

l-[e^- S(E)J g(Z-S) 

which is a result which vje have aready obtained in diagrammatic 
and mean field approaches. 

These three are the principal vievjpoints, basically 
equivalent, though the formal mathematics appears different. 

The final results therefore, not surprisingly, are the same 
in all cases. The reason v;hy we review all three in detail, 
is that we wish to examine an avenue of generalising the 
single-site CPA's to include cluster effects without sacrificing 
the Herglotz property essential for physically meaningful 
results. 



29 


2.4 Failure of ICPA 


The ICPA provides a reasonably good description of an 
alloy over a x-jide range of parameters. It X'/orlis x« 7 ell x> 7 hen 
the difference in orbital energies of the tX'/o components of 
a binary alloy, is small and becomes inadequate wlien 
.. LWL_ 




is large (Here is number of nearest neighbours) . 


Since ICPA is a single-site approximation it cannot 
take correctly into account the effects due to tx-jo or more 
sites. Firstly, it cannot correctly accoiunt for the randomness 
in the hopping-integrals involving two or more sites. In 
binary alloys the shape of the electronic density of states 
function as well as the band-widths depend on relative values 
of various interactions betv/een neighbouring atoms. These 
numbers can be verq;^ different for the two constituents of 
an alloy even though their atomic potentials are similar. 
Secondly, short-range order due to chemical clustering cannot 
be taken into account. Finally the effect of statistical 
clustering which is particularly important in the impurity 
bands, is totally absent. The impurity band shox-^s some fine 
structure (Payton and Visschsr 1967, Dean 1972, Alben et al 
1975) x< 7 -hich is lost in ICPA. This is a serious drawback 
because many physical quantities like conductivity, electronic 
specific heat, paramagnetic susceptibility etc. may be greatly 
affected by the structure in the impurity band (Vander Rest 
et al 1975) . 



The ICPA fails to give correct band- widths, even in 

the absence of off-diagonal disorder, the ICPA band width is 

considerably smaller than the true band-width. It is pertinent 

to quote here the ^ union theorem'' which states that if 

(Ep, /Ep ) and ,E^, ) are the lov;er and uooer band- edges 

1 2 ^2 

of the two constituents A and B of the alloy, then the band- 

edges of any alloy are (E. ,E„) where K = min (E„ ,E' ) 

• c i-c 1 2 1 

and E^ = max (E , E' . ) , From the theorem it follows directly 
®2 ®2 

that if E^ and E 2 are sets representing the energy spectra 
of A and B, then the spectrum of any alloy A , is 

"■ v— i. •• C 



The experiments on the magnetic properties of Cu-Ki 
alloys at large Cu-concentrations indicate that local-environment 
strongly affects the magnetic susceptibility of an individual Ni 
atom in the alloy (Kiimar and Joshi 1978) . In this regime the 
ICPA which ignores all environmental effect, does not ^give 
satisfactory results. The failure of ICPA for large W can be 
understood from the fact that for large VJ the electron mean-free 
path is small and in order to obtain a good estimate of the 
configurationally averaged Green's function we must consider 
scattering due to pairs, triplets etc. in a self-consistent 
manner. This is completely ignored in ICPA. 

The general, conclusion is that ICPA gives a good 
qualitative agreement with experimental data, certainly far- 
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superior to other earlier approximations* But in cases v/here 
local environmental effects are important/ we must tahe into 
account the cluster effects/ for a good quantitative agreeraent. 

2 »5 Generalisations of 1C?A 

From time to time suggestions have been put fc3rward 
about ways of generalising the ICPA to take into account 
the effect of clusters/ off-diagonal disorder and SRO. 

Extensive work has been done in this direction (Aiyer et al 
1969/ Seven 1967/ Capek 1971/ Cyrot-Lackmann and Ducastelle 
1971/ Freed and Cohen 1971/ Nickel and Krumimhansal 1971, 
Schwartz and Siggial972, Srivastava and Joshi 1973). 

Attemipts to generalise ICPA/ in the diagrammatic 
approaches have failed because of the enormous number of 
terms or diagrams to be summed. In the diagramm.etic approach 
if we retain skeleton diagrams which involve only the sites 
within a cluster of size n and reject all skeleton digrams 
V7hich involve sites outside the cluster/ we generate a 
cluster-CPA (CCPA) of a cluster of n-sites. Leath (l972) and 
Nickel and ^^ruramhansal (1971a/ 1971b) have given the procedure 
for pair-CPA (2CPA) . In the Fig. 2.3 v/e have shOTmi the 
cumbersome diagrammatic summation reproduced from the paper of 
Leath (-1972). The conplexity of the procedure is apparent. 

Any slight error in inclusion of any class of diagrams 
yields approximations that are non-Herglotz (Nickel and 
Krummhansal 1971a/ 1971b) . Any contem.plation of applying 



AA 



Fig. 2-8a 
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Fig. 2-8b 

Kapeat of fig. (2*8a) for the case of those diagrams (coltaas 1} tihlch would 
aps>ear in tim t matrix for scattering hy an Isolated pair of defects. 














Fig.2'8b 

Ibspaat of fig. (2*8a) for the case of those diagrams (coluaa 1} vhieh nouid 
appear in the t matrix for scattering 1^ an isolated pair of defects. 
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such a diagrammatic summation to bigger clusters in three 
-dimensional lattices/ appears intractable. 


In the .multiple scattering approach the system is 

partitioned into clusters/ each containing n- sites (Tsukda 

1972) and clusters are labelled as C. ,C„ ... C. ... just 

i z j 

as individual sites are labelled as R. / . . . R . . . . . 

12 j 

Let the vector (<R !,.,<R i) be denoted by {C^ 1 

nj+1 ‘ nj+n ' j . i 

and diagonal matrix with 


2 • , 1 • • • e . , 

nj+1 nj+n 


down the diagonals/ 


by e^ then 


V = I .J Pj = I V . 

J 


(2.29) 


where Pj projects the j-th element in the n-cluster 
sub space. 


In the multiple- scattering picture, initally the 
electronic wavefunction evolves according to an unperturbed 
Hamiltonian The scattering potential H “ Y ” — 

then switched on. If we a.ssume that S is cluster diagonal 
then electron suffers multiple scattering from clusters and 
scattering events from different clusters are represented. 

Tne T matrix relates the emergent electronic wave function, 
modified due to scattering, to the unperturbed incident 
wavefunction. as before (see equations 2.22 and 2.24 ) 

T = 7 t . and t . = T . (I + G. I t ) 

5 - 11 . 

Gj is an (nxn) matrix and . 

Ij = (Vj-Ej) [I - (Yj-E.) 
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The CCPA equation then reads 

=0 (2.30) 

The self-consistency condition (2.30) may be written 
interms of the matrix elements of the configurationally 
averaged Green' s- function, if we observe that the Green's 
function for the system v/ith the averaged medium everywhere, 
except the cluster (or cell) j can be written as; 

Gj = <G> + <G> 2.j (2.31) 

Therefore the condition (2.30) becomes 

<Gj> = <G> (2.32) 

which is similar to the ICPA condition. 

This equation when applied to simple model of random , 
binary alloy on a linear chain, with diagonal -disorder only, 
yields satisfactory results. The band-edges are produced 
more accurately than in ICPA. Ducastelle has shown that for 
CCPA with cluster-diagonal self-energy the approximate Green's 
function retains the Herglotz property of the exact Green's 
function which is a very favourable point, for the generali- 
sation in this form. But v/hen we apply CCPA equations to 
realistic situations the calculations become infeasible 
(Tsukda 1972). For this reason several other simplifications 
have been suggested. 



2.5*1 Central— Site App r oximat.ion 

Equation (2.30) and (2.32) are (nxn) matrix equations 

2 ■ 

and require solutions of n - non-linear simultaneous equations 
2 

in n - unkno’t'm matrix elements of the cluster-diagonal self- 
energy. To simplify the calculation of CCPA people have 
approximated the form of cluster diagonal self-energy. If we 
assume S to be site-diagonal v/ithin the cluster there will 
be only, one adjustable parameter in the effective Hamiltonian 
making it impossible to satisfy all the n -relations (2.30) 
or (2.32). VJithin these limitations the first approximation 
kno^'m as 'self-consistent central-site approximation' ( SCCSA) 
was due to Butler (1972) . In this approximation equation 
(2.32) is replaced by a single equation, 

<R i<G.>lR^> =<P. i<G>tS> (2.33) 

o — J o o „ o 

where is the orbital on the site at the centre of the 

cluster Equation (2.33) is equivalent to choosing a E in 
such a way that the averaged density of states per site at 

the centre of the cluster is consistent with the external 

0 

medium. This prescription though entirely acjiioc gives 
good results in the I'jeak scattering limit (Butler 1972, 

Brouers et al 1973a, 1973b) . VJhen applied to strong 
scattering region this approxim.ation begins to fail. At 
some energy region w’ithin tje bandG(Z) loses its Herglotz 
property. The density of states becomes many-valued and 
the sum- rule for the integrated density of states is violated 
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(nickel and Butler 1973) . These difficulties arise for 
simple cubic (Butler 1973) and diamond (Kumar and Joshi 1975) 
lattices for a 50—50 alloy with 6=1 where 6 is the 
modified relative separation in the well depths. In order 
to remove these difficulties Butler (1973) suggested a self- 
consistent boundary- site approximation (SCBSA) which is 
discussed in the next section. 

2.5.2 Boundary-Site Approximation 

In this approximation the consistency is required between 
the averaged diagonal element of ■ the Green's function at the 
boundairy site 

the cluster) and a diagonal element of the Green's function 
for the averaged medium i.e. 

<R^t<G^.>lR^> = <1^1<G>1R^> (2.34) 

For a linear chain Butler (1973) has shown that the SCBSA 
is exactly equivalent to CCPA. This is a surprising result# 
since the SCBSA also, like SCCSA# employs a scalar sel:f- 
energy. But this equivalence does not carry over to two or 
three dimensions. Kumar and Joshi (1975) have applied SCBSA 
to a 50-50 alloy with 5 = 1, on a diamond lattice. It does 
not seem to show any analytic trouble. In a detailed SCBSA 
calculation Kumar and Joshi (1978) have found that in the 
extreme split band case, some of the partial density of 
states do become negative in a small energy region.. Thus 
yielding Tonphysical results. 


G^.>iR^> (where b is a boundary- sate in 
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So far/ v/e have discussed the situations where only 
diagonal-disorder v/as present- Several attempts have been 
made to include randomness in the hopping- integrals (Blackmann 
et al 1971/ Shiba 1971, Poo et al 1973, Niizeki 1975, Kumar 
et al 1975). Here it is assumed that the diagonal term of 
the Hamiltonian can have either of the two values e^ and ep, 
and hopping integrals may have values ^aB 

depending on the configuration of the nearest-neighbour pairs- 
Even this model fails to produce more than simply an S-wave 
scattering ( Elliot at al 197 a 


2.6 Cluster- Embedding 


In addition to various approaches for cluster generali- 
sations discussed in the previous section, people have tried 
to embed large size clusters in scane non- self-con si stent 
effective (or average) medium e.g. immersion of a large size 
cluster in ICPa. or VGA medium (Diehl and Death 1979a/ 

Choudhry 1981). Here once the effective medium (ICPA or VGA) 
is obtained, one can write an exact Hamiltonian for a culster 
of n-sites as 


ntn) 
— c 


hCh) 

-ef f 


+ h 


(2.35) 


where H^^l is the Hamiltonian corresponding to a n-cluster 
-eff 

effective medium and h = "■ !?eff difference 

Hamiltonian finite only over the cluster. 
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where 

—erf 


The Green'' s function corresponding to 
= (ZI - )“^ : ■ 

= (ZI - 4^5, - 

= h - si?f y' dff 

= (ZI - 


( 2 . 36 ) 


each Orf ^ — eff (nxn) square matrix. 

The inversion of (I - G^£.p h) is that of a finite matrix, 
because h has non^ zero elements only in the clus'ter subspace. 


Here H 

— c 


(n) 


includes exact details of the interactions over the 
cluster while it sees the rest of the medium effectively. 
Consequently^ by choosing clusters of larger sizes, details 
of environment implanted can be gradually increased. If we 
choose a central site and its neighbours# the density of 
states at the central site obtained from the corresponding 
number on the diagonal of G^ is closer to the true 


Hamiltonian than the effective medium (VGA or ICPA) . 


For each possible configuration over a cluster, we can 
determine the density of states at the central site and also 
the configurationally averaged _Green' s function. Particular 
features in the density of state curve can be attributed to 
specific clusters. 


But this hind of cluster embedding Often gives rise to 
two kinds of difficulties: (1) the band-width remains incorrect 
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and (2) because the density of states still integrates to 
unity (one electron per site), this has the effect of 
distorting the structure specially near the band-edges 
to preseirve normalisation. 'Thus spurious structure may be 
obtained where there should be none. 

So, from the discussion of the previous few sections , 
it is obvious that none of the proposed cluster generalisations 
are satisfactory. We therefore, turn to an altogether 
different viewpoint of configuration-averaging, namely the 
augmented space formalism (Mookerjee 1973, Mookerjee and 
Bishop 1974, Haydock and Mookerjee 1974, Mookerjee 1975a, 
1975b, Kaplan and Gray 1976-1977, Kaplan et al 1980, Kumar, 
Mookerjee and Srivastava 1982). Basing ourselves on this 
forriialism, we shall develop an unambiguous, tractable generali- 
sation of the ICPA which retains the Herglotz property at all 
energies for all ranges of disorder (Kumar, Mookerjee and 
Srivastava 1982)* Off-diagonal disorder and short-range order 
shall also be incorporated. 

Recently in a series of papers(Kaplan and Gray (1976,77) 
and Diehl and Death (1979a, 1979b)) the authors have examined 
augmented space formalism and have generated self-consistent 
results, with off-diagonal disorder on linear chains in 
particular (Kaplan et al 1930) . There is a close and exact 
relationship of their algebraic equations to the graphical 
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method that has been adopted in our work. The graphical 
method has the added advantage that it is ideally suited 
in conjunction with the 'Recursion' technique (Appendix A) 
of Haydock et al (1972/ 1975) for the application to 
realistic models on three-dimensional lattices. The 
equivalent algebraic equations soon become tediously unwieldy 
as soon as the cluster size increases/ particularly on a 
three-dimensional lattice. This is the principal contri- 
bution of this v;ork and an indication of its future 


usefulness 



CHAPTER III 


AUGMENTED- SPACE PORI-IALISM CLUSTER COHERENT-POTEiTITAL 

APPROXIMATION 

3-1 Augmented- Space Method 

In this section we shall introduce a different approach 
to the problem of configuration averagingt-namely, the augmented- 
space formalism, first introduced by Mookerjee (1973) and 
subsequently developed in various aspects by other workers 
(see Chapter II). We shall describe here the construction 
of the augmented space and shall define the notation to be 
used throughout the remainder of the thesis. 

We start with the assumption that the lattice sites are 
randomly occupied by the atoms of the type A,B,C etc. The 
random variable e^^ then takes on the values e^, e^ or Cq 
etc. corresponding to the occupation of the site i. A 
set of values taken up by the random variables {e^} is 
called a configuration. A certain probability density 
P({e^}) is associated v/ith various different configurations. 
Assuming that different e^j are statistically independent, 

P({e.}) = n p (e.) (3.1) 

^ i 

where p^(e^) ig the probability density of the individual 
variables. This assumption ignores short-range order due 
to chemical clustering effect, which invariably leads to 



The more general 


statistical correlation between the {e.}. 

1 

case of dependent variables has been considered, by Kaplan and 
Gray (1976-1977). 

Haydock (1972) noted the fact that the probability density 
^i^^i^ satisfies the properties 

00 

p^(ei) >0 and / p^(e,) de^ = 1 

— OO 

These properties are specific to the imaginary part of a 
Herglotz function e.g. the resolvent of a self-adjoint 
operator in a certain Hilbert space. Taking clue from this 
Mookerjee (1973) suggested that a Hilbert space and 

a basis may be introduced such that p^(e^) corres- 

ponds to the imaginary part of the resolvent of a suitably 
chosen operator in tj)^, . If t f is a specially chosen 

member . of the orthonormal basis {lf(l)>l in then, 

n 1 

Pi(ei) = - I Im <f^| (el - M^)“^ ifQ> (3.2) 

where e e.+ lo 

1 

A prescription for construction of a suitable M for 
a given p(e) is as follows: 

For a given p(e), if one could find a convergent 
continued fraction of the kind. 




0 — » 
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Then a representation of the operator M is the tridiagonai 
matrix with a^z a 2 * etc. along the diagonal- and 

/ b 2 / • • « / etc. along the off-diagonal positions i.e., 

i a^ b^ 0 0 0 O 0 

! 

' b^ a2 b2 0 0 0 0 

I ^ 

I 0 ' b2 a^ fo2 . 0 0 O . 


in some basis { !f^>}. Of course, x\re have to restrict our^ 
selves to probability distributions all of whose moments are 
finite, otherwise a convergent continued fraction is not 
possible. Most of the physically valid probability distri- 
butions are of this kind, the only exception being the 
Lorentzian. The procedure involved in Haydock' s (1972) 
method of recursion for determining matrix elements of Green'' s 
operator is the inverse of this. 


For a binary alloy case with diagonal disorder only, 
e.g. ® represents the diagonal element of the 

Hamiltonian in a tight binding (or atomic orbital) basis. 

The diagonal terms {e^} form, a set of random variables 
described by 

e . = e'^N . + e^ ( 1 -K . ) 

where {N^^} are a set of randcxn occupation variables for 
the solute A in the solvent B’. In the absence of short- range 
order the various N^'s are independent random variables 
with probability densities 



4:3 


p(N^) = C6(N^-1) -i- (l~c) 6(N^) 


(3.3) 


I im ( = 


+ } 


n '♦O 


^ Im -r { -— ■ — + -- -— ■ —- i v?here N'=N.+in 


n-O* N^-1 


N' 


n ^ 


n -»0 N*-^- 


c(l-c) 
N^-(l-c) 


-} 


H I'" + !■ 
n -^o 


•V. , 

1 


N -a. 


where = c, = (1-c) and = c(l-c) . A representation 


of the matrix M is thus a 2x2 matrix 


I c ’^c(l-c) 

i’^c(l-c) (1-c) 

j 


(3.4) 


with eigenvalues 0 and 1. Its eigenstates I and 1^2^ 
describe occupancies corresponding to these eigenvalues. 


/Je now introduce the ^ disorder- field' ' . To each 

variable we have associated a vector space and a 

representation of a self adjoint operator corresponding 

to the bim.odal probability density of that variable in some 

basis If we had taken the set of eigenvectors 

{ |h^>} of M. as our basis, then each ih^> in (jj. would 

correspond to one particular value taken by the variable e^. 

In general any element in is a linear combination of the 

{ ih%}. (t>. with its basis and operator M. thus contains 

' n 1 
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all possible information about the configurations of the 

variable e; . 

1 

The product space § = H . can now be constructed 

i ^ 

v/hicn contains all possible states of the set {e^^} i.e^ all 
configurations of the disordered system. The basis in the 
product space is 

|f> = {lfH} {|f?>rix> ... (3.5) 

-1 

and if g (Z) = (ZI - M^) “ 

then P({e^}) =-| Im <f^iG(Z)if^> 


where | f = | f^> 0 i f 0 | f 0 ... 

(H ) (M„) (M.) 

and G(z) = g (Z) 0 g (Z) 0 g (Z) (xj • * 


(3-6) 


We call § the disorder- field. It is clear from the 
way it is constructed that it contains complete information 
about the statistical behaviour of the system. The infor- 
mation about probability of configurations is described by the 
operators defined on the field. 

We now describe the configuration averaging procedure. 

Let us first consider a function f(e) of one variable only. 
Its average is 

■CO 

f = / f(e) p(e)de 

— 'OO 

= / f(e) {~ I Im 

-OO Z-e+iO^ 

where p(e) = - Im (e+iO ) 



AssTJming novj that 


f(Z) is a function of a complex 


variable Z and has no singularities on the real axis, 

in the neighbourhood of the branch cut of the function 
M 


1 

2 Hi 


i 


M 


0f(Z) g (Z) dZ 

' CO 


(3.7) 


The contour over v/hich the integration is carried out 

M 

is taken around the branch cut of g (Z) along the real 
axis and not including any singularities of f(Z) * Since 
M is a Hermitian operator on ^ , we can -.-^rite 

OO 

f = - ^ <bf(s)dz <f^| / (z-h)~^dp(h) |fQ> 

—OO 

where we have written ^ (Z) = J' (Z-h) d^Ch) , £{h) being 

— OO 

the spectral projection operator of M. Thus 


f = f (h) dn(h) I fQ> 

— C50 


= <fo' ’^o^ 


(3.8) 


v/here f(M) is the same function of M as f is of e. 

The average has been expressed as a representation of a 
suitably constructed operator. 

The generalisation to a function of several variables 
is straight forward with one associated with each 

variable e^. In case f(e) is itself a matrix element 
of the kind say <r^ iP(e) i rj>, F eX state vectors 
if ,r.> = 1 ^ 0 ^ 0 *^i^ defined on the augmented product 



space ¥ - 3^ (X) n $ and the average becomes 

-A r “ ^ ^ 

F = <f^;rj_i[/ {Im G (h+iO ) } 0F(h) dh] lf„;r.> 

Vfe shall now consider the problem of configuration 
averaging of the Greenes function. For a substitutionally 
disordered, random, binairy alloy the Hamiltonian 

(with nearest neighbour interaction only) in the tight- 
binding basis can be v/ritten as 


H 


= I 


e.P. 

1— i 


ij 


V. .T. . 
IJ-IJ 


For a predominant diagonal disorder, we may assume 
and write 


H = Ze.P.+vlT.., HeM (3.9) 

- 1 1-1 ij-lj - 

where V is the weighted average of the overlap- integra 1 . 

This is valid for cases v/hen the bandwidths of the consti- 
tuents are roughly similar. 

Referring to the Hamiltonian of equation (3.9) with 
diagonal disorder only, the Green' s function can be written as 
G^^(Z,{e^}) = <r^ I (ZI-H)“^ irj_> 

The poles of G. . are situated at the eigenvalues of H. 

(M) 

These are distinct from those of g (Z) which are at the 
values that e^ may take. Using the averaging procedure, 
therefore, for every single random variable e^, we have 





(3.10) 


where |f^;r^> - | f^> (x) |r^>/ j f^> is given by equaricn 

(3.6) and 

S = - G. + V I T. . f>h I (3,11) 

—1 —1 , _ V,. _ V , X/ 

1 ij 

where 

i i * * • • • Q?) 51 j_ QS) i. 0^/ * ■ • 

(3.12) 

The expanded hamiltonian H contains complete information 
about the system and is defined on the expanded Hilbert ' 

(X) I 

space ¥ = J-| § where § = H . . Description of 

i ^ ; 

the quantum behaviour of the system is in the subspace ]-{ , 
while that of the statistical behaviour is in §, the disorder i 
space. The configuration averaged Green's function is a | 

particular representation in this expanded space W . The 
problem has thus been reduced to the determination of the ' 

resolvent of the Hamiltonian H in the augmented space ; 

Notice that up to this stage of configuration averaging no 
approximations are involved. Subsequently/ in the determination! 
of the resolvent in the augmented space/ herglotzicity retaining 
approximations may now be introduced systematically. In the ! 

next section we shall describe a method which will enable us 


to do just this 



3.2 The Graphical Technique 


The graphical method was introduced by Haydock (l972) 
for calculation of the resolvent. It v/ill be used to generate 
the cluster-CPA formalism for a general alloy system. 

By a graph we shall mean a set of ^ vertices^ connected 
by ■'links'. Let {ii>} be any complete linearly independent 
basis in which the Hamiltonian is described. The overlap 
matrix 

Sij 

is a unit matrix in case the basis is orthonormal. Matrix 
elements of an operator M are given by 

= <i|M|j> 

A, graph is associated with every basis {|i>}. To every 

element |i> of the basis set { li>} is associated a 

vertex v^, and to each distinct pair of elements (|i>y lj>) 

a link or bond . which may be directional (i.e. . 

3-3 "^3 

may not be equal to or otherwise. To any operator M 

there corresponds a graph consisting of vertices and those 
links for which 0. If x-je want to invert an operator 

M, then the contribution of each vertex v^^, is defined as 

k(v.) = 1/M., and that of the link *■. . as kCS-. .) +M . . 

A path of length N is defined as a series of 

ordered vertices and links i-e. 
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Q„ = { V . S, . . V . . = . v„ } 

IM 1 IJ J N' 

Contribution of a path is defined as 


k(Q ) = 
n 


k(v,) n k(£,,) 


vertices 


links 


ij 


A self-avoiding path of length N, Qj^={ v^ ^ v^ . . .v^} is a 

path such that v v with the only exceotion that v. 

■t' n m 1 

may be the same as Vj^z in which case it is a self- 


avoiding closed polygonal path. 


The renormalised contribution of a vertex Vj^ in a 


path is defined as 


k(v ) = I I k(Q^) 
N =0 


(3.13) 


where are all polygonal self-avoiding paths from v. 


n 


and back, on the graph from which the vertices 
have been removed. If M = G = then the 

representation of G in the same basis can be expressed 
■ (Hay dock 1972, Mookerjee 1979) as follows; 


OO 


G., = I I 

N=0 Q., 


(3.14) 


where are all self avoiding paths from v^ to v^. 


A use of graphical representation will be made, when 
approximations like CCPA are generated. 
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3.3 ICPA in the Augmented Space 


The Hamiltonian h defined on the extended space f 
= 3^ can be written in the site- field representation 


inf> as 


where 


H 


nf /mf ' 


M IS 

ff" 


6 + V 

~ffi nm — nm 
the operator in 


6 (3.15) 

ff" 

$ associated with the 


probability distribution and V is the ooerator in 3H . 

— nm 

The disorder field space states are specified by the set 

of states { |f^>} occupied in all individual fields 

the exact form depending on the probability distribution 

p(e). In general/ one chooses the basis in 4> ^ such that 

the representation of is tridiagonal so that p(e) is 

given by a continued fraction. For a birnodal density the 

/ for each n is given by equation (3.4) and its graph 
f f ' 

is a single link chain like 



v/here a^ = (1-c)/ a 2 = C/ b = /cil-c) . 

PInally the configuration averaged Green s function is the 
resolvent corresponding to the Hamiltonian H . 

Until now all equations are exact. Approximations are 
made while employing the graphical technique. The exact 



graph on which the calculation from equation (3.13) is 
formidable, is replaced by a simpler graph on which the path 
counting is( simpler and) tractable. 

Choosing for, M in equation (3.13) as (ZI-H) the confi- 
gurationally averaged Greenes functions G ( Z) can be 

nm 

determined by considering self avoiding non-intersecting 
paths in the complete augmented space §. Walks 

in the augmented space imply the following. An electron at 
a site labelled by n and field state i f > = (j?) 

i f (x) ...|f^>... can be induced by H either to make 

spatial hops to one of the near neighbours of n with the 
matrix element V, keeping the field state same, or it can 
remain on the spatial site, while the field at the site n 
changes according to the field at all other spatial sites 

remaining same. For | f > = 1 ^o^ ' ^o^ ... , 

G (Z) = <nf 1 (ZI-H)~^ imf> (3.16) 

nm 

In this sense the problem has been reduced to that of an 
ordered system but at the ejipense of a much amplified set 
of eigenfunctions. 

Starting from the vertex (Of) the first steps possible 
(fig. 3.1) are to the near neighboixrs in the direct 

space (shOTym by single lines) or to a different configuration 
state of the site 0 labelled by (0^) by a disordered field 
hop (shown by double lines) • 





/ 




«2'\ ”1/ 

\ / 

'' / 


w. 


(0,fQ) 




(Of) 


Fig. (3.1) 

if is the contribution of all self- avoiding 

non- intersecting paths from vertex 0 and back in the space 
14 , then the resolvent in an ordered system would be 

In three dimensions this will include the contribution 

J v^ ^ii^ where i are various neighbours of '0' and 
i 

corresponds to the resolvent calculated from a sub- 
graph in- which the vertex * 0 " has been removed, besides 
contributions of closed self-avoiding loops in 3-f . In 
the extended space, the configuration averaged Greene’s 
function G(Z) can be written as 

G(Z) = ^ (3.18) 

[Z-R(Z) - T(Z)] 

R(Z) is the contribution of all self-avoiding paths from 
the vertex 'Of^ and back that are either (a) entirely in 
the spatial part of the augmented space or (b) self— avoiaing 
paths in the augmented space which include field— hops but 



do not form closed loops. T(Z) is the contribution of all 
closed self— avoiding loops from ^Of^ and back in the full 
augmented space, connecting different spatial sites via 
field hops. 


The essential approximation is now made: in the augmented 

space f , there exist, closed self-avoiding loops which arise 

not because of any closed loops in the spatial lattice, but 

because of the field introduced. These paths contribute 

to T( and our approximation involves delinking all such 

paths such that T(Z) = 0. Once we have done that our graph 

consists of essentially that in the ordered system together 

(n) 

with an extra field hop with the off-diagonal element of M 
in the tri diagonal representation (b = / c( 1-c) in the case 
of binary alloys) as the link function. If starting from a 
vertex 'nf' v/e hop to a vertex then the subgraph x<i;ith 

the vertex 'nf^ missing is an infinite graph which is exactly 
similar to the original one. This property did not exist 
in the undelinked graph. A direct consequence of this 
property is 


= r, 

0fQ,0fQ ^0f,0f 


So 


(3.19) 


This will not hold in the absence of the type of 
delinking adopted above. 


Now we shall show that the above delinking procedure 
(Bishop and Mookerjee 1974) is exactly the 1 CPA. 



Let G (Z) be the resolvent, corresponding to the 
graph in the augmented space after delinking and let 

M. . = W. 5 . 6. 
ij 1 J/i+1 

be the tridiagonal representation of the operator M in 


each subspace of §» Then, 


G°^(Z) = 


[ Z-R^ ( Z-hf G^ ( Z) ) -wj G^ ( Z) ] 


oo 


* — -'ll '^-1 

o 11 

1 


P_. [z-VJ? G, (Z)l 


where 


•n 


( 3 . 20 ) 


= 5 2 . 


( 3 , 21 ) 


In the conventional ICPA formalism a self-energy S (Z) 
is defined by 


G^^(Z) = 


[Z-^(Z)-R^(Z-Z)] 

Poo [z-JU)] 


( 3 . 22 ) 


where g'"^(Z) is the ICPA- Green's function. The self- 
energy I (Z) is determined self consistently from the 
IGPA equation. 


OO y V 

J x-e(z) 

-oo l-G(Z) (X-5:(Z) ) 


p(X)dX = o 


( 3 . 23 ) 


so that 
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oo 


J 

— oo 


p(X) dx 
A(Z)-.X 


(3.24) 


where A(Z) 


1+G^^(Z) S 

gCP 


Nov 7 just as in equation (3,7) 
can be written as 


G^^(2) 


(3,25) 

from equation (3.24) 


.CP 


(Z) = <f^! [a(Z)I - M ]"^ If > 


o 


(3.26) 


Graphically therefore 
G^^(Z) = 


[a(Z) - wj G^^(Z)] 


(3.27) 


where 


gJ^(z) = 


[A{2) - 


(3.28) 


But A(Z) from equations (3.26) and (3.27) is the same 
DL 

as G (Z) from equations (3.20) and (3.21). This equivalence 
of two approximations has also been shown by Mookerjee (1974) 
using multiple- scattering diagrams. It establishes a direct 
connection with the delinking approximation of the diagrammatic 
approach to the CPA. 


3.4 Cluster CPA(CCPA) in the Augmented Space Method 

The essential approximation involved in the ICPA 
described in the last section is to delink all closed, self- 
avoiding paths involving both spatial and field hops in the 
augmented space. In the cluster - CPA, those self- avoiding, 
non intersecting paths are included which involve disorder 



field hops over the n- sites of the cluster. A self-energy 

renormalising the interaction between each pair of sites in 

tbis cluster can then be defined. The corresponding n x n 

2 

matrix for ^ will have less than n independent components/ 
which can be determined in a self-consistent way by an 
iterative procedure. 

The starting Hamiltonian is a tight- binding type: 
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+ Z I 

i/^j n 
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V. 

in/ jm 


T. . , 

— lb/ jm ' 


H e (3.29) 


where P and T are projection operators and transfer 

operators in ; n,m are band or orbital indices. Both 

the diagonal terms Ce. } and the off-diagonal terms CV. , } 
^ in iP/jm 

form a set of random variables described by 


e. = ^ N . + e? (1-N .) 
in in 1 in i 

V. . = V^N.N. + V®®(l-N.)(l-N.) + V^^{N.(1-N )+N.(l-N.) }, 

in/jm. nm 1 j nm i 3 nm i j J i 

(3.30) 

where {Nt} are a set of random ^ occupation' ' variables 
for the solute A in the solvent B. In the absence of 
short-range order the various N^'s are independent random 
variables with probability density 

p(hh) = c 6(N^-1) + (l-c) 6(N^) (3.31) 

Using (3.30) we may rewrite (3.29) in a more useful form 





H = + I I e N.P. + y T Y ( 1 ) 

“"B . ^ ni— in L L L 17 


. ^ n i-in ^ ^ V" 'N.N..T\ 

^ ^ 1.P3 n m nm 1 j — in^jm 


+ I l I v^^Un.+n.)t. 

i^j n m ^ J 

(3.32) 

A B 

where is the Hamiltonian of the pure solvent/ e = e — e 

““•D - n n n 


V 


( 1 ) ^ ,,AA , ^ ^B ,,(2) „AB „BB 


nm 


V + V - 2 V**"" and V'‘“' = V''“ - V' 
nm nm nm. nm nm nm 


The Hamiltonian is noxf in a form appropriate for the 
augmented space formalism. Corresponding to each there 

is a operator in a two-dimensional vector space 

A representation of in a basis { ivh/ ivh} is 


jC c(l-c) } 


1/2 


{c(l-c) } 


1-c 


1 / 2 ' 


( Hooker jee 1973) , 


v/here c is the concentration of the solute A. The augmented 
soace Harriiltonian is 


H = H 0 I + I e 

n 1 


I P 0 0 I 


(i) 


+ yy I T. . (x) 0 

nm . _2 .—in/ jm '■-..2 — — 


nm 


i/j 


+ y 1 T. 

L i/j 


+ (X! 


s "H 0 ^ (3.33) 

I ' in the above equation indicate identity operators 

f kV 

in all sub-spaces ^ except those superscripted. 


(i) _(ij) . 
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The configurationally averaged Green fiinction is then 
given by 

(3.34) 

with 1 n ^ 1 ^ 0 ^* Note that to this stage the 

relations are all exact. 

The density of states then follows directly: 

<n(E)> = -(1/7T) I Im <G. . (E+in)> ^ 4 - • (3.35) 

“ in/xn n-* 0 ‘ 

In the language of Kaplan et al . (l98C)/ it is the off- 
diagonal terms of which are responsible for the 

creation of pseudofermions or configuration fluctuations 
above the augmented space ^'ground state'' 1 ^; the 

presence of off-diagonal disorder pseudoferraions are simul— 
taneously created at two sites r^^ and r^ by the term 
fx) . 

Short-range chemical ordering may be described in an 
elementary way by a Markovian pair type— between neighbouring 

-St ^ -s 

sites r. and r. 

1 J 

p(Nj;lS!^=l) = {l-a(l-c)} 6 (N_.- 1 ) + a(l-c) 5(Nj) 

p(N.;N.=0) = ac 6 (N.-l) + (1- ac) 6 (N.) (3.36) 

3 Jl J J 

The single parameter a then is a measure of short-range 
ordering. If a = 1 then there is no short-range order and 

the alloy ds purely random# If 1 < cx < (l*"c) 


then . AB 
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type of clustering is favoured. The maximum value corresoonds 
to a situation when the alloy is not really random at all/ but 
forms a superlattice of A and B constituents. If a < 1/ 
then AA or BB type of clustering is favoured. 

This type of description of short-range order may be 
rather crude. Here we are dealing with a ^'quenched''' short- 
range order, that is, a short-range order that is not governed 
by thermodynamics and is consequently temperature independent. 
In more realistic situations one has to go through a more 
sophisticated formulation involving concentration-concentration 
fluctuations . 

Having formally performed the configuration averaging 
exactly, we shall now proceed to generate the approximation 
methodically, keeping in mind the important Herglotz property 
which must not be violated at any stage. We shall adopt a 
graphical formulation. It is well known that calculation of 
Green functions on connected graphs, using the nested path 
technique always produces Herglotz results (Haydock 1972, 
thesis). This methodology has been described in great detail 
by Bishop and Mookerjee (1974) and Mookerjee ( 1975a, b) * 

The graphs may be generated by the method described in 
section. (3.2). The essential approximation now involves 
ignoring or delinking the contribution of possible self- 
avoiding paths. The relation of these polygonal paths to 
the CPA has been discussed by Bishop and Mookerjee (1974) and 
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MoolcG£'jG 0 (19'75]d) h3.s dsinonstjrstsd th©ijr r'slHtion to ttis 

scattering diagrams* The philosophy of approxi— 
mation is as follows. 

We consider exactly all paths that involve only the 
spatial part of the augmented space, i.e. . Such paths 
are connected to the lattice structure and must not be 
approximated or delinked without destroying some of the 
feature of the londerlying lattice. In three and two dimensions 
the recursion method helps us to deal with this infinite set 
of paths. 

-3 jn - 4 

If we want a CCPA with (r^/r 2 / • • • within the cluster, 
then all polygonal paths that do not completely lie in the 
subspace 3 ^ and involve these vertices and none others must 
also be exactly accounted for. These are related to the 
multiple scattering within the cluster. 

Larger polygonal paths v/hich involve sites within and 
without the cluster but which do not completely lie in the 
subspace ZH must be delinked. The type of delinking depends 
very much on the environment we wish to immerse our cluster in. 
Immersion in a VGA medium yields the travelling pseudofermion 
approximation (generalization of the ATA) of Diehl et al 
(1979a). If we wish to immerse in a ICPA medium, we simply 
delink all these larger polygonal paths into a Cayley tree 
(Bishop and Mookerjee 1974) . The result is not self consistent 
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and in particular the band-width is the same as that of the 

ICPA. If we wish to embed our cluster in a self-consistent 

cluster medium^ we partition the lattice in Tsuhada type 

clusters/ delink the environment itself into cluster paths 

and self— consistently calculate the effect. We shall describe 

this in detail in the next section. It is instructive to 

note .at this point that this graphical technique is exactly 
.'§iV)h<3>v- 

equAv a J ^. B .:fe to the partitioning idea of Kaplan et al. (1980). 

3.4.1 The Self-Consistent Medium 

The calculation of the self-consistent medium in the CCPA 
is of central interest and is the most important contribution 
of this work. To illustrate the procedure let us first 
consider a 2CPA. We shall indicate the generalisation of 
bigger clusters subsequently. 

Note first that after the delinking procedure according 
to the prescription described in the previous section/ the 
resulting graph has the folloxving structure; each bond in the 
original lattice has/ in augmented space/ an ^'octagon'' 
decorating it (Fig. 3.2a). The eight vertices of the 
octagon corresponds to the eight different configurations of 
the two— site cluster (2 x 2*^) . Because of the cluster 
delinking/ no octagon belonging to a bond has any connection 
to another/ corresponding to another bond/ through an augmented 
space link. This exactly implements the idea that within this 



Vf = 2CV^^ic^V^’l Vqb 
V2=V^2)^C (1-C)V^^VYb8 
V3 = 2 {1 'O^fk (1 - Qf v<^VVbb 
t C1-C)[v%-C) 

rrPc)[v^^lcv^’^3 

C (1-C) 

yC(t-C)e 

5«b*’C« 

= ®B+(t-C)e 

Fig.3-2a 



represents an octagonal * =0^ 

decoration as in FTg.3-2a 

Fig. 3-2c 
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corirelated. scattering from three or more 
sites belonging to different bonds is ignored. The final 
erfect is that each bond is renormalized as shown in Fig. 3.2b. 
This gives rise to a matrix self-energy 



L. 


The aim is to calculate E_ self-consistently. It fully 
describes the self-consistent medium in the 2CPA. 

The bond renormalization procedure follovjs in two steps: 

:(a) The decorating octagon is not an isolated - octagon 
involving only the two spatial sites (0 and 1) . If we look 
at the vertex A(i.e* If^ in the notation of Mookerjee 1973), 
the bond AE (If^ to Of^) belonging to the octagon is only 
one of the Z bonds emanating from A to the Z nearest 
neighbours of 1 . The other remaining Z-1 bonds to 2f^,3f^,... 
also hang on to the site A. Moreover these bonds themselves 
in turn have their own octagons decorating them (Figi3.2c)- 
Thus the bond AE is itself immersed in a self-consistent medium. 
The same is tme for the bonds PG and HK of the octagon. We 
may take this into account by saying that the medium renorma- 
lizes the bonds AE, PG and HK. Once the renormalization is 
accounted for, the octagon is effectively isolated. The 
renormalized vertices Oq{Z) and links are foiind as 

Let us divide the lattice as follows into two 


follows. 
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subspaces (l) an unrenormalized bond AS, which v;e shall call 
system 1. It has a Hamiltonian v^(T ~ + T ). 

(2) A lattice L, which is the original lattice minus the 
link AS, and in which all bonds and sites are renoirmalized 
by • This has a Hamiltonian 


S*' ^ ^ P. + I. (z) Y T. (3) A linking 

Hamiltonian = E^(z) I + T^. + + T^^) . 

' ^ -P 

If we now use the notation X ^ denote the inverse 
of the operator X in the subspace y 

= £j_G = [ZI - fjCint) -j- P^ ^ 

where = [zi - ] ~-2 

with Pj_ + £2 = !• (3.37) 


Examination of the eqn. 3.3 7 immediately indicates 

that the effect of the rest of the lattice hanging on to the 

( 1 ) ( 1 ) 

bond AE is to change the Hamiltonian H to H + £ where 
the self-energy £ is given by 


a 


AA 


— '^eE ~ I 


k,me .K 


Ton 


’ae - °ea - H 4^,’ 


( 3 . 38 ) 


denotes a- nearest neighbour of A, 


(b) Let us now work on the renormalized, isolated octagon. 
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As before, let us divide the octagon into two subsystans 
a bond B (Of, if ) , the rest of the octagon M« 


H . , = II H. . 

ijeOf,lf ~ i3yof,if 

H = z Z (H. . + H. .) . 
iSB jSM 


As before, we have 

2® = £b e Eb = - S® - S®’ i 5“]'"® 


This immediately shows how the octagon renormalizes the bond B. 
We have 


= 

oo 




^BM ^ -m 

% ^jk ^0 


~ ^01 


M rSM 


jVO k;.^! 


(3,39) 


Equations (3,38) and (3.39) together provide a self-consistent 
set of equations for the calculations of the self-energy 
The various Green functions involved can be easily calculated 
on a three-dimensional lattice using the recursion technique 
of Haydock et al. In the case when we are using energy 
dependent effective Hamiltonians, the recursion has to be done 
separately for each value of energy* 

The cluster generalization to arbitrary large clusters 
now follows exactly the same philosophy. Equations (3.38) and 



(3.39) remain the same. The only modification is the 
augmented space unit. For example^, in the 3CPA the augmented 
space -unit is the 24-gon (Pig. 3. 2d) (Mookerjee 1973) and 



This has four independent components/ but note that both 
(3.38) and (3.39) for this case also yield four equations; 
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(3^40) 


we shall be interested in the diamond lattice^ where the 
nearest neighbour cluster consists of 5 sites* Here the 
decorating lattice is a 160-gon (Mookerjee 1974) and the 
self energy is still four independent members 
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As the cluster size increases/ the augmented space graph 
itself becomes large and the augmented space unit which 
renormalizes the cluster also becomes large. For clusters 
of Size N/ this unit is a (N x 2^) — gon. However/ the 

number of independent members in still remains relatively 
small/ and the can be calculated by the recursion method 

with great facility even on quite large augmented space units. 
This fortunate marriage of the above method with the recursion 
technique reduces the calculation of quite large CCPAs to a 
tractable problem . 

The self-energy S_ / hence the CCPA is cluster diagonal 
a property essential to preserve Herglotzicity (Butler 1973) 
The effective Hamiltonian is invariant with respect to cluster 
translations- However the exact effective Hamiltonian is 
translationally invariant. In contrast in the CCPA' s the 
self— energy at the cluster— centre and cluster— periphery are 
indeed different. Unless it is possible to construct clusters 
in which all the sites are equivalent (as in the 2CPA) this 
distinction is unavoidable and inherent in any cluster- 
approximation (Tsuhada 1972) . 

The self-consistent equations (3. 33) and (3.39) involve 
an -iterative solution for s_. We' could/ for example start 
with s°(z) = S^p^/ i:°(z) = V-in / use (3.38) to generate 

a, then (3.39) to get and iterate the procedure. This 

method will be useful only if this procedure converges. The 
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calculations were carried out on diamond and b.c.c. lattices 
with one electron per site and off-diagonal disorder. The 
procedure converged rapidly/ except near band edges, 
where the convergence is not rapid but nevertheless no more 
than a hundred iterations viere required in any case. The 
results are discussed in the following section. 

Results and Discussion 

(a) Diamond Lattice 

In this section x^?e shall present results for the 
density of states for a single s-orbital per site model on 
a diamond lattice. In order to compare our generalisarion 
to the self-consistent cluster C?A with previous ICPA 
calculations, we first present these ICPA results in 
fig. 3.3. 

Figure 3.3a shows the ICPA density of states per site 
in units of Ti~^ (-imaginary part of Green's function) for a 
dilute 90-10 alloy with small and large separations of the 
v/ell depths of the constituents. As expected, the latter 
case shows an impurity band. The characteristic feature 
of the ICPA is the lack of structure, particularly in the 
impurity band. The two— peaked structure in the majority 
band is characteristic of the diamond lattice and not a 
feature of disorder. Computer experiments, however, indicate 
peaks in the impurity band associated with clusters. Figure 
3.3b presents ICPA results for a concentrated 50-50 alloy 




Energy 


Fig. 3-3a 
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(ii) 

c =0-5 
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Vv'ith sm3.ll 3n<3. Isrgs sspsirstions of tlis wsll oeDths* Ag 3 in 
ths ICPA exhibits 3 l3c]<; of stouctuires other th 3 n those 
3rising from the ordered diamond lattice. 'The Van Hove 
peaks of the ordered lattice are rounded due to disorder. 
These four ICPA results will form a background against which 
we shall compare our cluster-CPA results^ to ascertain the 
effects of clustering and off-diagonal disorder. 

As mentioned in the previous section/ the self-energy 
self-consistent equations are practically useful only if 
the itera'bive procedure converges. Fortunately/ for all the 
subsequent cases not only did the procedure converge but/ 
except for the edges of impurity bands/ the convergence was 
rapid. In the exceptional cases the number of iterations 
required was large but tractable. For the sake of illustra- 
tion/ the self-consistent Zq(E) for the alloy 

later illustrated in fig. 3.6a are presented in fig. 3.4. 

The effect of clustering is apparent in X^(E) which has 
considerable structure compared with the smooth E^(E) 
characteristic of the ICPA. The effect of off-diagonal 
disorder and clustering is strikingly displayed in its 
effect on I.,(E). In the ICPA !.(£) is a real constant 
and is the average , value of the overlap parameter 
V - c^V^^+(l-c)^Vgg+ 2c(l-c)V^g/ whereas in the cluster 
generalisation it is a complex f'unction of the energy ano 
has considerable structure. 
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The first set of results illustrated in fig. 3.5 
considers only diagonal disorder in order that we may 
examine the effects of statistical clustering alone’. 

Figure 3.5a shows the density of states per site in 

-1 

units of n for a concentrated 50-50 alloy with small 

separation ( = 0.5) . Only half the density of states is 
shown as it is symmetric around the band centre. We first 
note that in this 'weak'' scattering regime clusters do not 
introduce much structure. The bond and cluster CPAS are 
very similar. Some structure is seen in the cluster CPAS 
and is in good agreement with the results of Mookerjee (1975b) 
which were not self-consistent. Secondly, the density of 
states near the band centre actually increases as we go from 
ICPA to CCPA/ indicating that clusters inhibit band separation 
For ICPA *5 = 1 is sufficient to separate the two bands# while 
for the CCPA 6 >1. This is also indicated by the moment 

method of Desjonqueres and Cyrot-Lackmann (1973) and the 
self-consistent boundary site approximation of Kumar and 
Joshi (1978). 

Figure 3.5b shows a 50-50 alloy with large separation I 

(6 = 1) . There is now greater structure in the CCPA results. ; 
The characteristic peak due to states localised on isolated I 

A atoms and two satellita humps due to stares localised on 
bonding and antibonding states on aB pairs show up. These ^ 

results are in good agreement with earlier work of Mookerjct; 
(1975) and various boundaxy site approximations (Kumar et al 1975; 
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Kiomar' and Joslii 1978) • The self-consistent central— site 

of Brouers et al (l973) suffers from violation 
of herglotzicity at E = 2.25, leading to unphysical 
discontinuity of the density of states at that energy. There 
is no band separation in the CCPA. 

Figure 3.5c' shows results for the 50-50 alloy v/ith a 
still larger separation ( £ = 2). The features described 
above show up even more clearly. Further, even in the CCPA 
there is now an unambiguous separation of the bands leading 
to a gap at the band centre. 

Figure 3.5d shows a dilute alloy (90-10) with large 
separation (6=1). This is the 'strong scattering' regime. 

The majority band is not much different from the pure host 
with the twin rounded peaks of the diair.ond lattice. The 
impurity band, however, shows considerable structure compared 
with the ICPA. The three-peaked structure mentioned earlier 
is nox'j even more prominent. For these parameter regimes the 
boundary site approximations start to show violation of 
herglotzicity. For 6 = 2, for example, around E = 3.2 and 
between 3.46 and 3.48 branch cuts appear off the real axis- 
This immediately restricts the use of these approximations 
precisely in those recfim.es of strong scattering where cluscciring 
has considerable effect. However, because of the graphical 
methodology and recursion techniques used in our v/ork, herglo- 
tzicity is assured in our results for all energies and 
disorders. 






Pig» '3*5(3 also shows the effect of Markovian short-range 
order of the type described in the earlier section and chara- 
cterised by a parameter (X.cx < 1 enhances the probability of 
precipitation of A clusters in a B background. This is 
clearly exhibited by the enhancement of satellite hump due to 
isolated A clusters, a- >1 enhances AB type pairs with a 
consequent rise in the central peak structure due to such I 

pairs. 

The next set of results introduces off-diagonal disorder ; 
in the overlap parameters. For any realistic model of an 
alloy it is ultimately essential to introduce off-diagonal i 

disorder. It is interesting to note that most previous 
tractable CPAS involving off-diagonal disorder necessitated 

to be either the arithmetic or the geometric mean of i 

and Vgg. This is an unnecessary restriction and ignores the j 
effect of chemical affinites in bonding. Our formalism, on 
the other hand, has no such restriction.. To start with we I 

shall take to be the arithmetic mean of and Vgg, so 

that we may compare with existing results but subsequently 
we shall also present results when this is not so. 

Fig.3.6a shows dilute (90-10) alloy with a large separa- 
tion (6=1) and off-diagonal disorder Vgg = -1.0, = -2.0 

and V = -.1.5. These results are to be compared with fig.3»5d* 

AB . 

Note first that because of the larger values of and 

The majority band resembles that of the host 


the band is wider 
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The structure near the lov/er Van-Hovs peak seems to be a 
feeture of of f“*clJ.aQonal c3.isorc3.er enci is also oresent in 
cases when <5 is small .(fig, 3,8). The striking new feature 
is the structure in the impurity band. Not only is the 
impurity band wider (as a conse.quence of larger overlaps 
betv/een impurity and host and impurity and impurity) , but 
there is a structure in the impurity band. 

In order to analyse the origin of these structures/ and 
to make sure that these structures do not arise due to 
spurious numerical errors, we shall examine the partial 
densities of states of specific cluster configurations 
embedded in the CCPA medium, as well as the corresponding 
eigenenergies of these isolated clusters. These details 
will unambiguously identify the origin of the structures in 
the impurity band. . The exercise is of importance since we 
have no completely reliable self-consistent CCPA results in 
this regime to compare our results with, certainly none with 
off-diagonal disorder. 

The configuration of the cluster corresponding to a given 
partial density of states in the impurity band is displayea 
alongside in fig. 3.6b. It is evident from figure 3-ob(i) — 
(vi) that the band between E 2L 3.20 and E o.60 is due to 
clusters of either host surrounded with impurities or vit-e 
versa. Its width is associated with V^g. On phe ocher hand 
fig. 3.6b (vii)-(x) the band above E 21 5.60 and stretching up 
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to E 6.38 is essentially due to clustering of impurities. 
Its width is related to The isolated cluster eicfen— 

energies also confirm this viev/ point. This is a specific 
feature of off-diagonal disorder clustering and cannot be 
expected to be reproduced in any ICPA in which the distinct 
effects of and cannot show up. The small 

peak at E = 5.80 is seen as due to isolated impurities, 
while the bigger peak around E = 5.60 arises due to clusters 
in which the host atom is surrounded with impurities. These 
are to be compared with the central peak and satellite humps 
in the impurity band of fig. 3.5d (without off-diagonal 
disorder) . The peak in the upper impurity band is due to 
large clusters of impurities (fig. 3.6b (vii)-(x)). This 
peak is small because the probability of such clusters is 
also small. 

Two isolated cluster eigenvalues at E = 7.279 and 
E = 7.636 lie outside the CCPA impurity band (shown by broken 
vertical lines in fig (3.6a) ). These arise from the cluster 
of the type shown in fig. 3.6b (ix) — (x) . A self-consistent 
larger cluster theory would eventually include these in the 
tail. Such isolated structures were also seen by Kiamar and 
Joshi (1978). 

Figure 3.6c shows the effect of Markovian short-range 
order. As before a < 1 diminishes the probability of unlike 
pairs. As a result, for a = 0 . 75 , the impurity peak due to 
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isolated A surrounded with B diminishes considerably. On 
the other hand oc, > 1 increases the probability of precipi" 
tating clusters -with unlike atoms, thus increasing the 
corresponding impurity peak. Ihis effect shows up in all 
cases illustrated in these figures. 

The partial densities of states and the short-range 
order effects serve to identify the origin of the structures 
seen in the impurity band. 


Figure 3.7 shows the density of states of a dilute 90-10 

/ 

alloy with off-diagonal disorder, but with a much weaker 
bonding. The qualitative features of the alloy in fig. 3.6a 
are retained. However, it is illustrative to note that since 
the overlap is weak, the peak due to impurity surrounded 

by host is not broadened, but remains a 6 -function-like 
spike around E = 4.3. 


In all the foregoing cases both 2CPA and the 5CPA curves 
are indicated. As expected, in all cases the 5CPA shovrs niore 
structure than the 2CPA. 

Fig, 3.8a shows a dilute 90-10 alloy with small separation. 
Apart from the structure near the lower Van Hove peak tnere 
are no striking features- Neither is there much difference 
between the 2CPA and 5CPA, This is to be expected in this 
^weak scattering'' regime. 
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Figure 3.8b shows a concentrated 50-50 alloy with small 
separation. This is to be compared with fiq. 3.5a. Unlixe 
the case of purely diagonal disorder the density of states 
is not symmetric around the band centre. This is due to 
the different bandwidth of the host and impurity. 

Finally, in figure 3.9 we have the last situation of a 
concentrated alloy (50-50) with large separations and different 
off-diagonal disorders. There is no symmetry around the 
band centres and there is a structure in the impurity band. 

(b) B.C.C. Lattice 

In this section we shall present the results for 
the density of states for a single s-orbital per site model 
on a b.c.c. lattice. The aim is a subsequent study of b.c.c. 
sd-metallic alloys like g-brass. These b.c.c. metallic 
alloys provide six sd^ orbitals per site. Although, we 
have considered here only one orbital per site the generali- 
sation to several orbitals per site is straight forward. 

Figure 3.10 shows the ICPA, 2CPA and CCPA density of 
states per site in the units of tt for a dilute 90-10 a,Lloy 
with <5= 0.75. In ICPA there is no separate im.purity band 
but in 2CPA and CCPA we observe a separate two peaked 
impurity band. Also there is an increase in the band-widrh 
from ICPA to 2CPA and CCPA. 







Figure 3.11a shows the ICP^ and 20?;^ (without off- 
diagonal disorder) density of states for the dilute (c=C.l) 
alloy but with larger separation in well-depths (6=1-0). 

This is the strong scattering^' regime. Here even in ICPA 
v/e observe a separate impurity band. Figure 3.11b depicts 
the 2CPii (with off-diagonal disorder) and CCPA density of 
states for the same alloy. Here again there is an increase 
in the band-width as compared to ICPA. The effect on 
majority band is very small and a single peak structure of 
the majority band is characteristic of the b.c.c. lattice 
and not a feature of disorder. Both 2CPA and CCPA show a 
two peaked structure in the minority (impurity) band. In 
order to identity the structure in the impurity band, the 
partial densities of states are shovjn in fig. 3.11c and 
3. lid. The peak around E '^ 9.6 is identified V 7 ith an 

isolated impurity ii in a B background (fig. 3.11c) while the 
peak around E 8 is identified with the cluster where AA 
interaction is dominant (fig. S.lld) . 

Figiare 3.12 shows the effect of Markovian short range 
order of the type described in the earlier section and 
characterized by a parameter cc, a < 1 enhances the probabi- 
lity of precipitation of A clusters in a B- background. 

This is clearly exhibited by the enhancement of the l-st 
peak in fig. S.lld. a >• 1 enhances type pairs which is 
exhibited by the enhancement of the second peak in fig. 3.11 
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Figure 3.13 shows a concentrated 50-50 alloy with 
' ~ O.dU. inliKe the ICPii case of ■ovirely disconal disorcler 
here in 2CPii. and CCPA the density of states is not s'^cr.T.atric 
around the band centre. This is a feature of the off- ' 
diagonal disorder because of the different band-widths of 
the host and impurity. 

These preliminary calculations on a simple but non- 
trivial three-dimensional model establish first, that the 
augmented- space formalism in conjunction to the real-space 
recursion techniques is a practical and useful method of 
generating physically valid (Herglotz) self-consistent 
cluster generalisations for resonably realistic models of 
alloys including off-diagonal disorder and short-range order. 
Secondly, the results themselves are illustrative of the 
effects of clustering, off-diagonal disorder and. short-range 
order on the density of states. We have so far considered 
only one orbital per site. The generalisation to several 
orbitals per site is straightforward. This work provides 
us with a sufficient theoretical basis for application to 
the physically interesting situation of III-'/ temary alxoys 

like Ga In- ^t-s, Go/is Sb- etc. These also form diamond 
c 1-c c i-c 

lattices, but there are four sp^ hybridised orbitals per site. 
This work has been reported in the next chapter. 



CHAPTER IV 


APPLICATION OP THE CLUSTER-COHERENT'-POTENTIAL 

APPROXimTION TO III-V SEMICONDUCTING ALLOYS 

4.1 Introduction 

The preliminary results of calculation on a simple 
model for random binary alloys discussed in the previous 
chapter were encouraging and provided us V7ith sufficient 
theoretical basis to apply the method of self-consistent 
cluster-CPA to more realistic situations of III-V alloys. 

The III-V semiconductors GaAs, InAs,GaSb^ InSb/ GaP 
and InP form zinc-blende structures (fig* 4.1). A zinc- 
blende structure consists of two interpenetrating f.c.c. 
sublattices/ the anion (e.g. As) and the cation (e.g. Ga) 
sublattices/ which are displaced from each other by a 
vector T = (i i i)a, where 'a* is the lattice constant. 
Consequently each anion in the crystal is surrounded by a 
tetrahedral arrangement of four cations and vice-versa. A 
schematic picture of the flattened network’ of the stiructure 
is shown in fig. 4.2. 

Pairs of these semiconductors mix in all concentrations 
to fonn ternary substitutional disordered alloys. The 
sublattice of the common ion is '’ordered^ while that of the 
other ion is occupied completely randonly. 
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Zinc blende crystal structure 

Fig. 4 -1 


O ANION 
0 CATION 
a (0,0,0) 
h (l/4,l/4,l/4) 
t (1/2, 1/2,0) 
d (l/4,3/4,”l/4) 

® (0,1/2, -1 /2 ) 
f (-1/4, 1/4, -1/4) 
g (3/4,3/4,1/t) 
li (1/2, 1,1/2) 

(-1/?, 3/4, -3/4) 
(l/8,l/0,l/0) 

2 (-l/a,-l/S,l/8) 

{-1/8,1/3,-1/B) 
(l/8,-l/B,-l/S) 



Cordinataa ax« ia '' 
tha tinita of a, 
tha lattiea oonataa 
of tha f.e.c, lattti 


Connectivities and coordinates in a flattened 
diagram for a zinc blende structure 

Fig. 4-2 
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The aim of the vrark reported in this chapter is two fold. 
Firstly, the direct application of the formulation developed 
in the earlier chapters to a realistic system of interest. 

The basic Hamiltonian is set up in the bond-orbital basis 
with parameters fitted to experiment or in the pseudo-bond- 
orbital basis within the ab-initio chemical pseudo-potential 
formalism (Choudhry and Kookerjee 1980, Choudhry 1981). 
Secondly, as we adopt basically a tight binding view-point, 
we attempt with the help of the recursion— technique (see 
Appendix A) to develop a completely 'real^ space method of 
generating the various constituent Green-functions. 

In our former aim we have been successful, as reported 
below. In the latter, however, the attempt has been only 
partially successful. Of course, very accurate k-space 
methods exist for calculation of the Green functions - for 
example, the ray- integration technique (Chen 1977) , which 
was further extended by Choudhry (Thesis 1981) to deal with 
off-diagonal Green function elements too. Green functions 
with vacancies (as required in our method) may also be 
derived from these by isolated impurity techniques. However, 
our aim was (to quote Heine) — throwing out l^pace (solid 
State Physics, Vol. 35) . However a practical snag appeared, 
that to get results of accuracy comparable to that of k-space 
technique, one has to carry out the recursion at least upto 
twenty to twenty- five steps. To make this feasible it requires 
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use of real space symmetries and an efficient way of extra- 
polating the asymptotic part of the continued fraction 
coefficients* Both tnese vjere beyond the knox'^-how of our 
group hgre. If we compare a seven step continued fraction 
version of the density of states per bond (in the units of 
7T ) for GaAs and InAs (fig. 4.5a,b) with a ray integrated 
version of density of states per atom (Chen and Sher 1978) as 
shown in fig. 4.6 then it is found that the overall features 
are reproduced with a s-like part/ an internal gap, followed 
by a p“like part consisting of roughly two parts in a ratio 
of 1:2. However, the early trxincation of the continued 
fraction produces spurious oscillatory dips. Since, the 
completion of this work, it has come to our notice (Kelly, 
private communication and Solid State Physics, V'ol. 35) 
that efficient programmes have been developed that take into 
account real space symiTietries (Kelly, Solid State Physics, 
vol. 35, p. 244) and truncation schemes which claim to 
achieve accuracies equal to the k-space methods. The 
author proposes these as further extension of this work. 

However, he would like to emphasise the above, as a further 
justification, for starting to develop a fully real-space j 

technique, wedded to the ideas of the augmented space method 

j 

introduced in Chapter III. 

fi 

4.2 Hamiltonian [ 

me work by Harrison and his coworkers (1973-74) has i 

demonstrated that the aystematics of bonding of various | 
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physical properties can be described by bond orbital model 
(BOM) in such covalently bonded systems. In this model 
the basis functions are centred on bonds connecting each anion 
and the nearest cations. These basis functions are constru- 
cted as linear combinations of sp^ hybridised orbitals on 
adjacent atoms. This construction results in bonding and 
antibonding orbitals. Instead of atomic orbitals as basis, 
the bonding and antibonding orbitals form an equally valid 
basis. We may partition the Hamiltonian as H = H +H + H^ 
where and project on to the valence (or bonding) 

and conduction (or antibonding) subspaces respectively and 
is the hybridising term. 

-P 

Sv = - Sv- sJc 2, 3 '' 

Ge = - Sc ] ° 

The valence band desnity of states is given by 
- — Im I . within the bond orbital model, it is assumed 

7T —OCX, OCX 

LL 

that H is small, and so we may completely neglect the anti- 
—cv 

bonding part and work exclusively in the bonding subspace. 

In case there is a significant coupling between bonding and 
antibonding states# the bond orbital approximation is invalid^ 

and bonding basis is incomplete to describe even the valence 
band only. 

We use the chemical bonding orbitals as our basis set 
and using the notations of Chen and Sher (1978) we fix the 
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anion f.c.c. sublattice as the reference lattice on which 
each lattice point is represented by a lattice vector j. 

The four bonding orbitals surrounding it may then be 
represented by ija>, a = 1,2, 3, 4 centred at points specified 

•Mb 

by j + where =(l,l,l)a/8 (-l,-l,l)a/8. 


T3 = (-1, 1,-1) a/8 and x^ = (i,~i,-i)a/8 

v;here a is the 

lattice constant. The bond orbitals themselves are composed 

of symmetric linear combinations of the hybridised atomic 

orbitals |a> and i c> centred at the anion and cation atomic 

1 - 

sites respectively i.e. | ja> = [_!&> + ic> J = iac> where 

n 

is the normalization constant ’^ 2 ( 1 + 3 ), S being <aic> 
the orbital overlap of anion-cation pair, and 


ia> = I |s^> + *^1 |p^> 

1 ^ 

|C> = ^ |s^> + ^ lp^> 


( 4 . 1 ) 


where |s > and ip > are the s and p orbitals on the 
0 . 0 

anions while is > and those on the cations# 

The p-orbitals are taken to be directed along the bond in 
the above expressions .When expanded in terms of a bond 
basis set { lta>} the 30 M Hamiltonian takes the form 

H = - ^ , I j,a> <j,a|H i j'£x'> <j'a' i 

ja,j a 

If we consider bond to bond interactions upto second 
nearest neighbour bonds, there are five distinct interactions 
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which are defined below (see fig. 4.2). 

D ; The bond energy <ja|H^l^a> 

3 . 

; The matrix elements between adjacent bonds with an 
anion in comm.on. 

Y ^ % The matrix elements between adjacent bonds v.dth a 
cation in comraon. 

; The matrix elements between parallel second nearest 
neighbour bonds- 

^2 5 The matrix elements between non-parallel second 
nearest neighbour bonds. 


The corresponding Bloch basis functions can now be 
constructed from the bonding orbitals in the usual manner: 


■A 1 

|k,a> = ~ 

/N 




(4.2) 


where N is the total nximber of anions or cations in the 
crystal and k the wave vector v/ithin the Brillouin zone. 


The secJ-lar equation can^ therefore, be written as 

(4.3) 


det iE(k)l “ (k) 1 — C 


where 


H^^'(k) = I uTL exp ik.CjtV 

3 oJ 


(4.4) 


has been obtained by using the transformation of tases 
defined in equation (4.2). Here, one of the site suffixes 
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on the Hamiltonian in the real space has been fixed as the 
origin . 

In the above description , represents/ for e>:amDle 

oo 

a 4 X 4 matrix representing the energies and interactions 
between the four bonds emanating from the anion site at the 
origin. Similarly by assigning different values to the 
lattice vector j/ interactions between other bonds can be 
brought in. 

Chen and Sher (1978) have shown that the inclusion of 
third-nearest neighbour interactions/ does not change the 
band structure significantly. SO/ retaining the interactions 
upto second nearest neighbour parallel and nonparallel bonds 
only, the distinct matrix elements are D/ T^/ Y;|^ and ^ 2 * 

In the work of Chen and Sher these form independent parameters. 
The k dependent 4x4 Hamiltonian matrix under this approxi~ 
mation is therefore, 

Hii = D+ 2 D 2 [cos( 2 x+ 2 y) + cos( 2 x+ 2 z) + cos (2z+2y) ] - 
= D+ 2 D 2 [cos( 2 x+ 2 y) + cos(2y-2z) + cos ( 2 x- 2 z) ]- 
H 33 = D+ 2 D 2 [cos( 2 x- 2 y) + cos( 2 y- 2 z) + cos ( 2 x+ 2 z) ] . 

= D+ 2 D 2 [cos( 2 x- 2 y) + cos( 2 y+ 2 z) + cos ( 2 x- 2 z) ] . 

H ^2 = 2{Y1S cos(x+y) t 2 y 2 cos(x-y) cos(2z)!- 2iYlA sin (x+y). 

H ^3 = 2{y 1S cosCx+z) + 2 Y 2 cos(x-z) cos(2y)}- 2iYlA sin (x+z) . 

H. , = 2{Y1s: cos(y+z) + 2 y 2 cos(y-z) cos(2x) H 2iYlA sin (ytz) . 



H 

H, 

H 


23 

= 2 {YlS* 

COS 

(y-z) 


2^2 

cos(2x) 

COS 

(y+z) } 


2iYlA 

sin (y-z) 

^24 

= 2{y1S 

cos 

(x-z) 

+ 

2^2 

cos(2y) 

cos 

(x+z) } 

-f" 

2iYlA 

sin 

(x-z) 

‘34 

= 2{y1S 

cos 

(x-y) 

-f 

2^2 

cos(2z) 

cos 

(x+y) ) 


2 i Y A 

sin 

(x-y) 


where ylS = 2 Y lA = (y^ - y-^)/2 and X/y^z 
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stand for j k^a/jk^a, - k^a- respectively . 


4.2.1 The Parameter Fitting Procedure 

Chen and Sher have adopted the empirical parameterization 
approach. ITiey have obtained the values of parameters by 
fitting energy separations at specific syminetry points on t-he 
Brillouin zone r /X/L etc. The results are given in the 
table 4.1. The parameter D has been determined emprically 
using the band energies and the experimental values of the 
photo-electric thresholds measured by Shevchik et al (1974). 

Table 4.1 




Parameters (Mev) 



Material 

D 

Y IS 

YlA 

°2 

^2 

GaAs 

-11.1687 

-1.5187 

-0.4250 

0.3406 

-C.1CB4 

GaSb 

- 9.6437 

-1.2687 

-0.3750 

0.3906 

-0.0594 

Gap 

-10.8125 

-1.3625 

-0.3250 

0.3812 

-0.1188 

InAs 

-10.0562 

- 1.3937 

-0.5125 

0.3844 

— 0 * 06 5' 6 

InSb 

- 9.5625 

-1.3500 

-0.3875 

0.2937 

-0.0563 

InP 

-10.2750 

-1.3125 

-0.4125 

0.2625 

-0.0375 
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In the table 4.1 we have shown only an optimised set of 
values. Often, variation in the values of a parameter 
obtained by fitting to experimental data, for a given material, 
is more than that over dif rerent materials. Chen and Shsr 

finally impose an energy dependent Lorent 2 ian broadening and 
vary the parameters in order that the broadened peaks thus 
obtained best resemble the experimental spectra obtained by 
Ley at el (1973). 

The procedure adopted by Chen and Sher is arbitrary and 
its validity is also not clear. Further the bonding orbitals 
are not orthogonal. The arguments of Chen and Sher, willing 
away the overlap integrals, are of doubtful validity. 

The alternative and. much more satisfying procedure is 
the abr- initio chemical pseudo-potential approach of Anderson 
(1969) and Bullett (1975), as performed on GaAs, InAs etc. 
by Choudhry (Thesis 1981). Within this formalism the bases 
used need not be orthogonal. It is sufficient that they 
fonTi a linearly- independent set. It must be mentioned that 
the chemical pseudo— potential calculations were not self- 
consistent /as such their accuracy is also in some doubt. 

However, we were not interested in this aspect of the 
problem and for our CCPA calculations, we have used, the 
parameters of Chen and Sher. 
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Table 4.2 


Bond lengths, bond energies, and work 
functions for Ili-v semiconductors 


Material 

s Bond length 

(A°) 

Bond energy 
(e.v.) 

Work- function 

(e.v.) 

GaAs 

2.45 

-5.6687 

5.5 

InAs 

2.61 

-4.7562 

5.3 

GaSb 

2.65 

-4.7437 

4.9 

InSb 

2.81 

-4.7625 

4.8 

Gap 

2 . 36 

-5.1125 

5.7 

InP 

2.54 

-4.5750 

5.7 

4.3 Rec 

ursion Calculation 

of Green . Functions 


The alloys under consideration are of the kind Ga^In,__^As 
as a common anion alloy and GaAs^Sb^_^ as an example of 3(a) 
common cation alloy. These alloy in all proportions and 
therefore c may be anything between 0 and 1. Each of the 
six materials of table (4.1) have nearly the same bond 
length (table 4.2). Hence we shall ass^lme that the intera- 
cation parameters corresponding to the pure materials are 
the same as those in the alloy. The vacuum state sets the 
zero of the energy scale* 

In an alloy the question of interactions between two 
different kinds of bonds arises* For example in the case of 



Ill 


alloys/ a configuration of the hind 


o h 

Ga. In 


ijia;^ be a part of the alloy, it is demonstrated below that 
the interaction i^i2 between bonds 1 and 2 is the average 
of s corresponding to pure Ga-As and pure In-As crystals 

as a reasonable approximation. 


Let the bonds 1 and 2 -be described by the q^lantum states 


|b^> and 

lb2> 

respectively. 


Then 

lb^> 

1 

” N 

( l^s^> + 

1 > ) 

and 

lb2> 

1 

“ N 

( | As2> + 

tln2» 


3 

in terms of sp hybridised orbitals at various sites directea 
along the respective bonds and N is the normalisation constant 
^<■^2, assuTTied to be the same for the GaAs and InAs bonds* 

We are concerned here with H ^2 i-e. 

<b^ = -2 +<Ga2.' ^ ^ 


{ 1AS2> + iln2>} 1 
= — „ [<As^ |HiAs 2^ + tGa^^ iHiAs2^ 

+ <As^lH|rn2> + <Ga^lHlIn2> ] 

The first term is vanishing because of the orghogonality 
of two different hybridised orbitals at the same site. 
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The middle two terms vroulcl be identical in the unailcved 
case due to the geometrical symmetry of the location of the 
o j- -Oi tci 1 '-j iii'/olve'J» Hence this cart of the exoressio'.' is 

ecTi.iivalent to the average of these corresponding to the 
pure cases. ihe last term is small conmared witli the 
middle two since it involves more distant orbitals. The 
approxiraation therefore involves the neglect of the variation 
in this term over the constituent miaterials. 

In order to compare the results of cluster-C?A forraalism 

with those of ICPA and VGA v/e restrict ourselves to the one 

soecific alloy i,e, Ga^ ^-As, However the procedure 

0.5 0.5 " 

for other III-V alloys is quite similar. 


In order to get the electronic density of states per bond 

->*•1 , „ , 

scaled by n " , we have used the recursion method of Haydock 

(1972 ). The parameters used are those as used by Chen and 


Sher. The VGA parameters are given by the average of those 
of the constituents weighted by their respecti'/’e concentrations. 
The VGA and the self-consistent ICPA results are shown in fig. 
4.7. The ICPA results have been obtained by using the ICPA 
equations 


c(e^- e “ 2) 


( e - e - S ) (1 -c ) 

= O, 


l“(e^ - e - S) g(z- S) l-(eg - e - S ) g{z- i) 


A 

where 


c e^ + (l “C ) eg. 


The VGA diagonal Green function Is substituted as a neroeth 
approximation and the quadratic equation in ^ thu^ ob^-ained 
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is solved. The choice of the correct solution is ensured 
by checking that X o as 6-0 and the imaginary part 
of a remains negative throughout the band. At each step 
new g(E-a) is calculated and substituted back in the abcve 
sguation. -Lhe process continues until sslf~consistent 
solution is obtained at a fixed energy. The energy is 
varied in regular steps across the whole of the valence band. 

A comparison of VGA and ICPA curves shows that the CPA 
rounds off the sharp peaks present in the VGA. 

4 *4 Cluster-*CPA formalism 

The single- site CPA completely ignores the effect of 
clusters and cannot adequately deal with off-diagonal 
disorder. Chen and Sher have treated diagonal disorder in 
ICPA while off-diagonal disorder in VGA. This is not a 
consistent procedure and the randomness in the off-diagonal 
elements/ particularly in the nearest neighbour overlaps/ 
is significant. Here an attempt has been made to improve 
upon this drawback. 

When GaAs and InAs are alloyed together they form a 
random alloy. One of the f.c.c. sublattices is ccmpletely 
occupied by the As atoms/ while the other f.c.c, sublattiee 
is randomly occupied by the Ga or In. The analysis therefore 
resembles that for the randan binary rather than the ternary 
If we go to the bond-orbital basis and label our 


alloys. 
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bonds by the positions/ then the bonds are either as— 'Ga 
or AS- In type randomly, ihe pseudo-Hamiltonian cicsely 
resembles that of an Anderson model. 


H 


1 


I O.P, + I I 

■ i/jSyV i/jeb/ -iJ 


^ t . T. . + Yo F , , ‘T. . 

i/j s i\f' ^ l,je 

where (or denotes nearest neighbour transfer 

between tv/o bonds with an anion (or cation) in common. 
^;N|'^(or /\/’^) denotes nearest neighbour bonds with a cation 
(or anion) in common. /V’ ' ' /\f ^ denote next nearest 

parallel and non-parallel bonds (see fig. 4,2) 
labels a bond |r,n>. 
are both random/ so that the model contains both diagonal 
and off-diagonal disorders. The second nearest neighbour 
parallel and non-parallel bond interactions D 2 and Y 2 
respectively are treated in VGA, denoted by D 2 and ^ 2 * 


Here 


The set {D.} and the set CV. .1 
1 1 j 


We may v^rite H interms 


of the occupation number H . for 


the (say) Ga-As bond. That is if I'b - 1, the bond is a 
Ga-As one and zero, if it is not. Then 


or 


D. 


D. 


^GaAs(jj^) + D^^^'®(l-Mj_) 


D'^^i + D®(l-N.) 


where D*"" and D® represent bond energies of Ga-As 


and 


In-As respectively 
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Now 

1.3 

and 

t;t( 3. ) 

V . J can be wrrtten 

as. 

( c ) 
V. I 
IJ 

11 

N. 

X 


• 


y ( a ) 

ij 

AiA 

= r, 

X 0 

M . N 

1 

j.:- 

+ 

-brid-ii 

la X J 1 j 

Here 

AxA f 
^ 

1 CL 

or 

BB ^ 

7^ ^ ; represents 

first 

nearest nsighdoour bond 


interaction between two Ga-As (or InAs} bends with As as 

common ion and j\ (or ) represents that betv/een two 

1 c ic ^ 

nearest neighbour Ga-As (or In-As) bonds with Ga (or In) 
as comnon ion. r" represents nearest neighbour bond 
interaction between Ga-As and In-As bonds with As as cormor- 
ion. so H can be written as 




+ I T. . -i-(r. 


K.., = 


V . -IJ 
IJ 


B3 ) y m 

- 2 ' Ij Aj 


n® y p 4- y r,n(o) y rj(s) 

D z ?i + z_ ‘ ^la ij 


IJ 


IJ 


-1- D 


....E3 


I T. . w r 

Lt - 5-1 


2 


2 « * 
^ ij 


is the pseudo-Hamiltonian for the pure InAs case; 


D,- D,.. and 

■ A . tb. 


AA ^BB _ „ AB 
3a “ "^la '' 3a la 


„AB ^BB 

3a = 3a ' 3a 


^ 1 /,„AA 

2C 2 


BB' 


= ~ '■3c Ic^ 


(4.6 ) 
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Here In-As has been used as sci /enb and Ga-As as 
solute. We micht "iust as v;ell reverse their roles. 


:arn 


rderinc due to cherrdcal 


inn a 


= co(n^-l) (i~c) 6(K.) 

lere c is the concentration of Qa-As. 


(4.10) 


VJe new go over to the augiTienred s'Oace formalism. This 
has been described in great detail in Chapter III. Using 
the same notations introduced earlier^ the augraented space 
Hamiltonian can be written as 


.B 


I 

i 


+ V, y T. . 

la r . —XI 
11 


+ (V^^+ v„„) 


XI I-. 


'ij 





(i) 


0 iClj) 

(i) 




I 

+ ( y^- Y®®) 1 




(j) 


-ij 


vjhere i' 


,-Ji) 


''c(I-g) I 


! 


(4.11) - 

(4.12) 


b^(l-c) 1-C 

L 


and ii eH 6c) 


The. con figurational Iv averaged Greenes function is 
given by 


<G^^(Z)> = <i, T°i (21 - H ) ^1 j, 


(4.15) 


o 6?) 

where ' ] t > = n 1 v 

i 


is the augmented space configuration 
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"ground state' (chapter III). The density of states is 



<n ( 

E)> = 

i 

Im, 

<G^. ( S-fin) 

• 

•I* 

A 




i 1 

o 

X Jl. 


'Tlie 

CC?,-. m 

.3.Y t)e 

■generated 


) using th e r - ^ fn o d cl c gp' 

des 

1d0O. 

in Ch 

apeer 

X X a « 




The c 

lalcul 

ation 

Or ^ 

lelf-consist 

:ent nvediuxn is of central 


interest. j.n craar tc rliustrate ■c'ns oroce t.irs '-"e shall 
first consider a pair-G?A (2CPA) . If we apply the delinking 
procedure described in Chapter III, the resulting graph has 
the following structure. 

Each two oonds attached to a common ion in the original 
lattice have in aucjmented space an ''octagcn'' clecoratina it. 
(fig. 4.3a). In the case of Ga^In^_^ As, there are two such 
octagons, one corresponding to two bonds with a cation in 
coramon and the other corresponding to two bonds with an anion 
in common. As a starting point, we treat the interaction 
between two bonds with a cation in corrm.on in VGA, and denote 
it by i.-e- there is no cation disorder and the intera- 

ction between tx-^ro bonds X'/ith anion in ccmuTion is taken into 
account exactly. We can do othervjise also, i-e. we can taka 
into account exactly the cation disorder while treat the anion 
disorder in VGA. Our choice is governed by ',-;hich disorder 
is more significant in our example. In order to take both 
types of disorder into account consistently, we must perform 
a 3CPA calculation. However in the present work v/e have 





perfonneci only a 2-;pA calculation cind a full 31PA and cchar 
cluster CPA calculations have been left for laser research 
work. The au-grienced space Hard 1 ton ian with hond disorder/ 
cue to a co'u.'.cn ar.icn onlp/ can be '.'risten as 
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(4.15) 


the eight vertices of the octagon (fig. 4.3a) correspond, 
to eight different ccnfig’urations of the tx*;o bond cluster 
(N X t ) . i/Je should notice here that each vertex of our 
graph nov; corresponds to a ■'bond' and each link to the off- 
diagonal Hamiltonian matrix element in a bond rather than 
■'site' basis (In fact it is the 'covering' lattice formed by 
a bond- vertex interchange) . 

Now- w^e apply the cluster delinking procedure described 
in Charter III. The result is that no octagon belonging to 
a pair of two bonds with an anion in comraon has any connection 
to another corresponding to another pair of bonds with an 
anion in common. This again implements the idea that within 
2CPA correlated scattering from three or more bonds is ignored. 
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The final effect is that each pair of bonds with an anicn in 
coTf/iCn is renorrralissd. This gives rise tc a laatrix -elf- 

il: I , " 

; o a ; 

i 

T = i > 

■ y y I 

Here i is the self-energy of t:ie bonds and r, , is the 

*-'■ i. 

overlap— self energy’' of a nair of bends '-/itli an anion in 
comnaon. rhe aina is to calculate self consistently* 

The renonaalisaticn cf pair of bonds with an anion in 
common follows in two steps: 

(a) The decorating octagon is not an isolated octagon 
involving only , the two spatial bonds (0 and 1). If we look 
at the vertex A (i.e. If^)/ the link AE(lf^ to Of^) belonging 
to the octagon is only one of the six links emanating from A 
to the six nearest neighbour bonds of C. Tlae other remaining 
5 bonds to also hang cn to the bond A. Cut 

of these 5 bonds/ 2 hang on to the anion and 3 hang on to 
the cation constitutincr the boxid 0 (see fig. 4 -3b). These 
bonds themselves intum have their ovm octagons decorating 
them (fici. 3.2c) . Thus the link AE is itself immersed in 
a self-consistent medium. The same is true for the links EG 
and KK of the octagon. We may take this into account by 
saying that the medium renormalises the link AE/ EG and HK. 
Once the renormalisation is accounted for the octagon is 
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effectively isolated. The renormalised vertices £^(t) 
and links CT^(Z) are found as folloi’S Let devide the 
lartice (of bor-fc/in^cc two sub spaces (l) An unrencrrr.alised 
link AS, which we snail call systera 1. It has a dar'ilconian 



(2) A lattice L, which is the original lattice niinus the 
link AE and in which all bonds and those links which have an 
anion in common are renormalised by Z . This has a 
Hamiltonian 


. rC 2 ) 


^ ^i J. . . 

jii ij^h./ . 


T. . + Y, 

— ij Ic 




-1 1 


+ Dj Z T + 7, I 


IJ 


2 e -ij 


and. a linking Hamiltonian 
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v/e nov 7 use the notation X ^ to denote the inverse 

of the operator X in the subspace y we get 
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Exarr.ination of the form of above' set of equations 
Jia'cely indicates that the effect of the resr of the 
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Here ^ denotes the neighbours (upto second nearest 
naignbour) of bond A* The different Green functions/ used 
on the above expressions are specified below; 

Go “ diagonal matrix element of Green function 
GlA matrix element of the Green function between two bonds 
with an anion in common 

GIC - matrix elemtent of the Green function between t%ro bonds 
with a cation in comraon 

G2P - matrix element of Green function between bfO second 
nearest beiahbour parallel bonds 
G2A - matrix elemient^of Green ^function between two second 
nearest neighbour non-parallel bonds* 

G3?A - matrix element of Green function between two third 
nearest-neighbour parallel bonds with their; closer 
ends adjacent to anions 

G3PC - matrix element of Green function between two 'third 

nearest-neighbour parallel bonds \^urh their closer ends 


adjacent to cations 
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Giir' — iTia"cri:'C sj-ement or Green function betv/een toro fourth 
nearest neighbour parallel bonds 
G4 a — ioiatrix element of Green function between two fourth 
nearest-neighbour non-parallel bends 
G5PA — matrix element of Green function betv/een two fifth 
nearest-neighbour parallel bends with their closer 
ends adjacent to anions- 

G5PC — matrix element of Green function beb'/een tv/o fifth 
nearest-neighbour parallel bonds with their closer 
ends adjacent to cations 

G5AA ~ matrix element of Green function between two fifth 
nearest-neighbour ncn-parallal bonds v/ith their 


closer ends adjacent to anions 
G5AC - raatrix element of Green function between two fifth 
nearest-neighbour non-parallel bonds i-fith their 
closer ends adjacent to cations 


(b) Let us no'w v/ork on the renormalised isolated octagon* 
As before let us divide the octagon into two subsystems 


a link 3 (Of^ If) and the rest of the acta ton M 


Then 
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This immediately shows how the octagon renormalises 
the link B. We have 
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equations (4.17) and (4.18) together provide a self-consistant 
set of equations for the calculation of the self energy _ 

The various Green's functions involved can be easily calcula- 
ted using the recursion technique. In case when we are 
using energy dependent Hamiltonians# the recursion has to be 
done separately for each value of energy. 

4.5 Results and Discussion 

Using the methodology described in the previous section# 

we have generated the valence-band density of states for 

Ga^ Tn^ _ As. 

0.5 0-5 



Fig. 4.4. displays the the electronic censir" cf stctac 
(IFGS) curve for Ga--4s with nearast-neight our herd inccracricr. 
caly, •■ 7 iuh eigne and nine step recursions- The teihiinazion 
of oha centinued fraction of coefficients is cased on 
Gershgorin thecrem (see appendix-A) . Hera uhere are too 
many cscillaticns in the EDCS-curve an'"' we can'not distinguish 
betv7ean tv*70 parts of the p— like part of the band- 

Figure 4.5a shows the SDOS-cum/es for Ga-iis with six 
and seven step recursions/ with second nearaso-neighbeur 
bond interactions also. Here we could net go beyond seven 
steps because of time and storage-space problems. Figure 
4<.5b shows the seven step recursion curve for InAs. 

If we compare the seven step recurion SDOS-curves for 
GaAs and InAs with the EDOS-cur^/es obtained by k- space 
methods (fig. 4. Ga, 4.6b) for these materials then we note 
that the recursion curves have almost all the qualitative 
features of the curves obtained by k-space methods: 

(i) There is one s-like part with a peak and a secondary 
peak simulating the k-space hump. This has 25 per cent 

of the total area ('under the full curve) under it- 

(ii) A central gap insensitive to hov; the continued fraction- 
expansion is teirminatad. 

(iii) A p-like part vjith one peak with 25 per cent of the 
total area under it/ separated by a dip 'from a section 
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tunately both of these were beyond the knovj-how of the group 
here/ w?ien the work on this prcblam was in progress- 

In fig. 4.7 vja have sho'-jn the 1C?A and VGA results for 

Ga^ ^Ah. It is to be noted that the sinole site CPA 

O.o O.D 

does not have much effect other than smocthening the VGA 
results to a non-significant degree. 

Figure 4.8 depicts the ICPiTL result . If we compare 
this with the ICPA result / then we observe that there is 
no change in the internal band gap. There is a slight 
increase in the band width and a slight enhancement of the 
s-tvpe peak/ which reflects the effect of InAs. The lov^er 
rj~tvpe ;peak gets lowered and broadened. 






The importan-c thing to note is that no new structure 
appears in 20?^^^. This is because of the fact that there 
is verr/ small difference between GaAs and InAs parameters- 
This is in accordance with our predictions based on model 
calculations on diamond lattice (Chapter III). Thera also, 
when the difference in 'parameters was small, no new structure 

In conclusion, we have presented a method of obtaining 
a self-consistent cluster CPA. In particular, we have 
applied it to ternary Ga^In^_^As and calculated the valence 
band density, of states incorporating off-diagonal disorder. 
This is the main contribution of this work. 



CHAPTER V 


CONCLUSION 

Oi.ir attempt to extend the theory of disordered alloys 
beyond ICPA has resulted into the development of a self- 
consistent cluster coherent-potential approximation (CCPA) , 
which can be applied to alloys of materials very different 
in mass# potentials and bond lengths. This is a significant 
development. 

Though, the problem has been formulated only for 
electrons# the same model and mathematical treatment carry 
over to different kinds of elementary excitations like 
phonons and magnons. , 

To sum up# the following are the achievements of this 

work; 

(1) A self-consistent CCPA based on augmented space 
formalism of Mookerjee, has been developed to take 
into account the effect of large clusters. 

(2) The analytical properties (Herglotzicity )of the Green 
function are maintained in ' the CCPA-theory at all 
energies# concentrations and degrees or disorder. 

(3) The CCPA theory takes into account exactly the effects 
of off-diagonal disorder and Markovian type short-range 
order. This is a significant advance over earlier work. 
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(4) llie CCPA theorv has been successfully applied to Iix— v 
ternary alloys which is the be gin in g of application of 
the theory to realistic alloys. 

However/ there are a few refinements still essential, 

over the calculations presented here. They are 

(1) The real- space methods for the calculation of different 
Green functions (G. .) need further refinement. In the 
calculation of density of states of pure GaAs and InAs 
we could go do\Am only upto seventh step in the 
recursion because of vdiich there are some spurious 
oscillations/ near the Von-Hove singularities in the 
density of state curves. These oscillations can be 
avoided by going down more steps in the recursion. 

This can be facilitated by the use of real space 
symrcietrias vhich will reduce the working tune enor- 
mously. Haydock (Solid State Physics/ Vol. 35) has shown 
that the entire basis l of tight-binding functions is 
unnecessary/ for the recursion method/ only those 
combinations which have the symmetry of the lattice and 
the starting state i 0} / (see appendix A) ccn tribute - 
This drastically reduces the 'effective'' basis and thus 
the time and storage-space required. 

Another approach is to take a finite cluster and 
use the 'super smooth' programaae of the Cambridge group 
to smoothen out the resulting histograms. 



Pertorming a CCPA calculation which takes into account 
the full off-diagonal disorder i.e» both anion and 
cation disorder in the case of III-V semi-conductors. 
This is possible again^ by the use of ^ super- smooth' pro 
grammes or by the exploitation of the symrri^tries of 
the underlying lattice- in the recursion programme. 

Instead of usincf pararnetrization procedure (which is 
very ambiguous) of Chen and Sher the chemical-pseudo- 
potential method may be used for the calculation of 
different nea.r neighbour interaction parameters in 
III-V semiconductors. 

The variation of bond length has also been neglected. 

In the pseudopotential scheme, the interactions can be 
recalculated at the interaromic distances corresponding 
to specific alloys. In this xvay, the approximation 
made for Ga-As-In kind of near neighbour bond intera- 
ctions can also be refined. This involves the 
calculation of somewhat more number of interaction 
parameters. 

The effect of conduction band on valance band has been 
ignored completely. Just as bonding orbitals yield 
Green functions for valence bands, anti— bonding 
orbitals determine the conduction band. The full 
band i.-a. both conduction and valence bands can be 


obtained by using a tomdc- orbitals as basis. 


The 



, tii-Coct oi conduction band, sacms to altcn the top of valance 
b'Hiii'l ''/hich is important for all transport properties. 

Though, we have concentrated mainly on the problem of 
determination of electronic density of states or <G(S)>, 
it is possible- to extend this foniialisin for the calculation 
of ditfer€ 2 nt electronic px'operties of the disordered systems. 

The formalism can be used to explain optical properties 
or the phenomenon of photo-emission in solids. Here the 
number of random variables will be more, but still the 
problem can be tackled and optical- transition probability 
can be calculated (Joshi and Mookerjee 1982). Also by the 
use of Ward^s identity and accompanying Bethe-Salpeter 
equation, vre can make predictions about conductivity and 
dielectric response of a disordered system. 

The effect of bigger, clusters can be included by per- 
form.ing a recursion calc^ilation in the full augmented space. 

So far we have confined ourselves only to crystalline 
systerrs v/ith compositional di.sorder. Amorphous alloys, 
where there is a departure from perfect crystallinity, but 
there is an underlying distorted network on which the atoms 
sit (Polk 197 1 , Connell and Temkin 19 7 4)» may also be 
talked by the use of 'distorted'', orbitals as basis. In 
such cases, Hamiltonian matrix elements will nave continuous 
rather than bimodal distributions, but the augmented space 
techniques are general enough to consider these too. 
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APPENDIX A 
RECURSION METHOD 


•me 'recursion' method developed hy Haydoch et al (1972, 
1075) facilitates the detemination of the Green function. 
Given the Hamiltonian of a system,, . in suitahle Hilhert^ space 
the recursion method sets up a 'basis in which the namritonran 
has a trldlagonal representation, from which the matrix 
elements of the Green function may be derived. Ecjuava en 
,lve„ the Hamiltonian matrix H, the recursion method generates 
, ,„itary transformation U such that 

U H u — Std 

1 r, H is tridiagonal, 

where Hmj^ 

,, H in a Hilbert 

het us start with a suitable Hamirtonxan H 

The ouantity of interest is the resolven 

Space on . 

also called the Green operator 


G(Z) = him ^ (ZI - H) 
Z-»E+iO 


-1 


. . 4- ctio:! mans the Hilbert space on to 

H is self-adjoint and maps -n ^ ^ 

. - “U or, is Hiu >. Now beginning 

14= 4 if la > Sv- 1 n 

Itaelf i.e. if I n , ,,, , and an arbitrary sequence 

with a oiven starting element 1 a„> 

/ la b , let us define a sequence of elements 

of numbers n^ 


{ 1 u > I by 

' n 


2 ' . _ ■ ' 'O' 

u > n-l n-2 ^ O 


iA-l ) 



Tlii;:; sequence terminates at n if lu .> is the 
. n“f*x 

riii.l.L vtK.’.tor (bj^=0) / otherwise v/e generate an infinite 
;u'.]i,KTic;'.,' of olcraents. The numbers are real if 

H is Herinitian. If the sequence of elements { is 


such that . <u^|u > = 0 for n in then equation (A-1) 

m . ^ 

immediately implies 


n 


, .2 _ 

<u |u > n <u , [u 

n ' n n-1 n-1 


<u lu > 
n ' n 

<u - lu .> 
n-1 ' n-1 


(A-2) 


If v/e take { |u^>'} as our basis of representation/ the 
Hamiltonian , H becomes tri diagonal v;ith the down the 

diagonal and .down the off-diagonal positions. 

If H is a non-Hermitian matrix (as is often the case 
with effective Hamiltonians arising out of the chemical 
pseudo-potential approach or any coherent potential, approxi- 
mation) then/ alternatively v/e choose an arbitrary set of 


complex numbers follows 


b^|u^ > = (H-a^) |u^> / c 


0 


'o'^l . 
for n > 0 

b lu > = (K-a ) |u > - c |u„_^> 
n n+1 — n n n n r 

Using equation (^3) it is easy to see that 


(A-3) 




n 


n I ' -n ' - n 'n-1 ' — ' "n 

2 

n ’ ' n' ■ n' • n ■ n ' 


< u„. 1 ( H-a^ f ( H-a^ ) ! >- l c. 


^ = 6 . 

if v/e assvime <u^|u.^> nm 


(A-4) 
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The reason why this method is 

Called recursion method 

now becomes elear. Starting from , . i, 

^ ^ Ham.il toman and starting 

element iu^> v/ith b^=0 we may , - 4 .^, -u t 

o o - ^aerate, with the help of 

equation (A-l)/ a^,b^ and thus determine 

in turn aj^,bj and thus |U 2 > and recursively. As a 

calculation technique this method 1^ economical. At 

each stage we need to store !u^> 

coefficients calculated upto that extraneous 

information is left out. 

Given it is simple to various matrix 

elements of the Green function, can evaluate the local 

matrix element 

G (E) = <01 CEI-H)~^ 10> 

OP . 

= < 0 1 ( jiI-K,pQ) i * 3 _f^ as we can/ 

we insist that 1 0> remains crs the f^j-st member of the new 
basis { i»s» ~ lG>. -^e off-diagonal element 

of the Green operator also be calculated in 

same way.' If we start with • , 


'V 


1 r 


[ |01 + !l> ] 


Here =<0i(EI-H) |1> and 


<%KSi-H) ^ !u^> - 2 [GoQhf. ] 


DO 


■O, 

= G, 


11 °ol = "^10 

-1 


GfoCE) = <Taol(EI-H) !Uq>>,q_(e) 


oo 


taut 
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Here the first term of the right hand . side of above 

equation can be calculated by continued fraction expansion 
1 

~ [ 1 0> + il>] and is the term which has 


with I u^> 


already been obtained by the same procedure with |u^>,= i0>. 
Hence we can calculate G^^(E) = (E) . 

X u ux 

Let D be the determinant of the matrix 

( n=0, 1 , 3, . . . ) be the determinant of the matrices derived 

n 

from (EI-H^^p after rows and column fraa 0 to n have been 
eliminated. Then if we wish to obtain Gqq(S) = 

\ 1 SG the Gauchy expansion for the determinant and write 


down a recursion relation for as 


D = (E-a )D ,,-b 


2 .. 


n 


n^ n+1 ‘"n n+2 


(A-5) 


Goo(E) becomes a continued fraction expansion from the 


repeated application of equation (A~5) 


D T 

G (E) = 2 = — — 

°° (E-a^)-b,DyD^ 


(A-6) 


E-a -br 
o 1 


E-a^-b2 


E-a2-b3 


b' 


n-1 


E-a -,-b„(: 


"n-1 n^D 


By 


n-1 

r the recursion method we can calculate iteratively 


(a, ,b ) from n=o to some N. 
b n 


The resolvent is then 



calculated by fixino 


in one of the ways described in 


N+1 


n 


N+2 


the? next subsection (a- 1) and inserting this value 
N— th level, of the continued fraction 


at the 


The density of states can be expressed ass 


n(E) 


■ Im 
Z"*E+iO 


4 * 


-l/n 


Z-a 



2 

1 





(A-7) 


The tridiagonalisation of H reduces the original 
system, to a fictitions linear chain of atoms with local 
orbitals with diagonal energies a^ and overlap elements 



(A-1) Termination of Continued Fraction 

9 

In most of the cases settle dOT-jn quickly 

2 

and tend towards a sym.ptotic values . By setting 

(a /b^) = (a for n > N i-je can evaluate the ram.ainder 

of the continued fraction expansion analytically. We write: 


G „(E) = 




E-a -b G_ 

CO oo OO 


^ %+2 ^ N +3 


and solve for G^,(E) yielding ■ 



must 


we reject the second solution of the auadratic as G (E) 

voiii.‘/li Us E -► oo ii'} Order that the resolvent preserve its 

analj/ bic character (of the form E ) there. It is clear that 

real . for real E outside the range ^^oo~^l^ool) ^rd 

(a^+2ib^l)., so that the continued fraction will be real-. 

Inside the same interval G_^(E) has a branch cut and we 

obtain a band of allowed energies. If we underestimate the 
2 , 

value of b^/ we will get spurious structure in the inteirval 

9 

(a -2|b l), a +2 ib i) but if v/e overestimate b the band 

CXJ CX3 QO OO CO' 

width is artificially large but the spectral v/eight in the 
spurious tails is extremely small (Haydock et al 1975) - 

2 

In the tables (A-1) and (^2) the values of (b.-,) for 
diamond and b.c.c. lattices (w^ith one orbital per site) 
I'espectively, are given. Here the well depths and overlaps 
are zero and -1 respectively. The well depth being zero all 
the are zero both in diamond and b.c.c. lattices. It 

is easy to see that in case of diamond lattice a^= 0 and 

9 2 
b = 4 v;hile in case of b.c.c. lattice a = 0 and b = 16. 

, Q0 CW WW 

The corresponding density of state curves are shown in 

figs- Al and A2 respectively. With the number of coefficients 

N = 18/ the density of state curves are fairly good. 
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Table A-1 

Recursion-coefficients for Diamond Lattice 


9 ,2 


IN 

1 

4 . OGO 

10 

3.723 

2 

3.000 

11 

4.210 

3 

5 .000 

12 

3.8 30 

4 

3.400 

13 

4.185 

5 

4.365 

14 

3.793 

6 ' ' 

3.747 

15 

4.216 

7 

4.292 

16 

3.817 

8 

3.663 

17 

4.148 

9 

4.352 

18 

3.873 


Table A-2 




Recursion-Coef f icients 

for B.C.C. 

Lattice 

— 2 

N 


N 

K 

N 

1 

8,000 

10 

16.47 

2 ■■ 

19.00 

11 

15 .58 

:3i," 

14.26 

12 

16.39 

4 

17.32 

13 

15.65 

15.01 

14 

16.33 


16 .8 3 

15 

15.70 

7 

O 

15.31 

16 

16.28 

16.61 

17 

15.74 

9 

15.48: 

18 

16.25 
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When we have several orbitals per site (e.g, four 
3 . . ' 

sp orbitals pe^r site in case of III-V- semi-conductors) 

2 ■ 

settle down so quickly/ as they 
incorporate the detailed information of intemai band - 
edges and singularities* Hence it becomes difficult to 
choose an appropriate ^ Ib^l }. This can be done in a 
reasonable manner by several methods (Solid State Physics 
Vol. 35 ) • We shall not go into details of these methods 
but shall describe a method based on Gershgorin theorem 
(Franklin 1968) which is, in general, suited for the 
estimates of (a^, !b^d) in case of isolated bands e.g. 
the valenoe bands of III-V semiconductors. 

The Gershgorin theorem states that for any Herrriitian 
matrix M, its eigenvalues E must lie within the union of 
teqions defined by the relations? 


xM - I IM. . d < E < M. . + I liM. .d (A--9) 

ii j ij - - -‘I ->1 


'ii ■ j ij 


for all i. Therefore by direct examination of any 

Hamiltonian matrix, v/e can determine upper bounds for the 

limits of the band, where the density of states is non- 

zero. Tills theorem gives very good estimates of (s^/ ^ 

when applied to the tri diagonal matrix (Kiplan 1978). We 

can thus determine the maximxim energy (E ) and the minimum 

f^nerov (E . ) limits for the band. Now knowing the fact 
^ min 


that 



v/e have 


1 

“ 2 

(E + 

E . ) 
min 

max 

1 

"" 16 

(E 

max 

E . ) 
min 


(A. 10) 


In case of valence bands of III-V semiconductors, the 
band width and the internal gap produced by a seven step 
continued fraction version, are in very good agreement, with 
those obtained by experimental and k- space methods (see 
Chapter IV). However, the early truncation of the continued 
fraction produces spurious oscillatory dips in order to 
preserve normalisation. This has been reported earlier by 
Haydock et al (1975) also. These spurious oscillations 
can be eliminated by going down more steps in the recursion, 
v/hich is possible by the incli „lon of symmetries of the 
covering lattice or by the use of finite clusters and 
subsequently using the ' supei?- smooth/ type of programme of 
the Cambridge- group (Solid State Phys. Vol. 35). These 
techniques have not yet been mastered by us and should be 
incorporated. 

The only thing that remains is the consideration of 
singularity structure of the density of states. In the 
complex Z-plane, the Green function exhibits all three types 
of singularities; simple poles, branch points and essential 
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Gingularities. The first few coefficients 

in the continued fraction expansion of Green function, are 

;.ru, tl'icient for the location of simple poles and the value 

of Green's function near singularities, but the position of 

band edges (branch points) and Von-Hove singularities 

(essential singularities) is determdned by the asymptotic 

behaviour of (a , ib i'') for large n. It is not possible 

n n 

2 

to tell the nature of singularity from first few (a^,ib^i ). 
Our termination of the continued fraction expansion after 
certain finite niimber of terms makes it impossible to 
determine the precise analytic structure of the singularities 
but the first few coefficients quite accurately distribute 
the spectral weight about the singularity (Haydock et al 
1975 ). 



APPEI^DIX B 


RESOLVENT OR THE GREEN'- OPERATOR 

The resolvent or the Green' operator G corresponding 
to the Hamiltonian H is defined to be 

G(Z) = (ZI-H)~^ (B-1) 

v/here I is the identity operator. It is defined for all 

values of the generalised complex energy variable, 

except on a domain viz. the spectrum of the Hamiltonian H. 

For a self-adjoint H, the spectrum of H lies entirely 

on the real axix in the complex Z-plane. For a finite 

system the resolvent has a finite number of poles (degenerate 

or otherwise) on the real axis. For an infinite system 

we can have (1) a set of poles (2) a branch cut i.e. an 

infinite, absolutely continuous (dense, compact) set of 

poles (3) essential singular! ties-lihe winding points 

(4) an infinite singularly continuous set of poles (dense, 

non- compact) as in the case of Anderson localised states in 

random Hamiltonians. The Green operator has poles at 

each spectral point , i.e. eigenvalue. Projection on to 

the corresponding eigen-state is the residue at that point. 

The diagonal member representation of G(Z) 

v/ith respect to any linearly independent basis ! ls>} is 

analytic everywhere in Im Z 0» Also Im ‘3gg(Z) < 0 if 

Im Z > 0 and Im G (Z) > 0 if Im Z < 0. If the spectrum 

ss 
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1 

of H is bound from above, clearly GggC^) 2 

,1 ixis In case H is not bounded from above, 
along the real axis. m 

,,P h:we to further assume this asymptotic behaviour of 

G (Z). Finally G(Z) , All these properties 

S3' 

define a "Herglotz" function. 

in some representation { , s» the Green function 

is given by. 


as 


CO 


G ^(Z) 
ss 


< S t C ZI-H) ^ 1 S> 


(E-2) 


= I 


<sin> <n}s> 

Z 

n “ 


, ■ eh-> complete set of eigen-states of H 

where { in» is the oomplet 

oorrospondlng to the eigenvalues (F,!, normalise 

^luy. Ihen the trace Of the oreen function matrix m 

this representation 


Y Y <^m> <nis> 


= ^ Z-E 
n n 

lateness Of the tesis set C,s>l is used. 

where the oompleten 

in case of tight-binding basis m a crystal,^ 

1 ^nts are identical and we can writ-e 
matrix elements 


gCz) = d 


ss 


(z) =1 I 

n 


E 


n 
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Since the eigenvalues are all real, the g(Z) 
defined along the real axis in the limit that the imaginary 

part of Z becomes very small i.e. 

1 


Liru g( E+in ) 
n ^0+ 


Lira 


n E-rin 


n 


* f ^ 

* I * 


_ i. y - iTl S (E~E ) ] 

- H ^ E-E^ 

E being real and (P represents the principal part 
so that 


- - Im g(E+in) =1- I 

” n-*0+ ■ " 


n 
n(E) 


where n(E) is the density of states per site. Ihis relation 
is very useful because taowlng n(E), the complete Green 

function at any value Z can be obtained from Kramer' s- 

Kronig relation 


1 


CO Im g( S + iO ) 


g(Z) = - n 


Z - E 


dE 


/ 

— OO 


n(E) 

Z-E 


dE. 


This relation is a oonseguanoe of the Kergolotsiclty 


of g( Z) . 



